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G> ■ We address nonperturbative power corrections to inclusive decay widths of heavy 

r-| ■ flavor hadrons in the context of the 't Hooft model (two-dimensional QCD at N^ — >■ 

^! cx)), with the emphasis on the 'spectator-dependent' effects, i.e. those sensitive to the 

OhI flavor of the spectator. The summation of exclusive widths is performed analytically 

using the 't Hooft equation. We show that the l/mq expansion of both the Weak 
Annihilation and Pauli Interference widths coincides with the OPE predictions, to 
the computed orders. Violation of local duality in the inclusive widths is quantifled, 
and the new example is identifled where the OPE prediction and the actual effect 
5^ I are completely saturated by a single flnal state. The qualitative aspects of quark 

hadronization emerging from the analysis in the 't Hooft model are discussed. 
Certain aspects of summation of spectator-independent hadronic weak decay widths 
are given in more detail, which were not spelled out previously. We also give some 
useful details of the 1/mQ expansion in the 't Hooft model. 
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1 Introduction 

The decays of heavy flavor hadrons Hq are shaped by nonperturbative strong interac- 
tion dynamics which, at first sight, completely obscures most of the properties of the 
underlying weak interactions self-manifest at the quark level. It is suffice to say that 
the actual hadrons, rather than quarks are observed in the final state. The actual 
dynamics of confinement in QCD to a large extent remains mysterious. Nevertheless, 
significant progress has been achieved in describing heavy flavor decays applying the 
formalism based on Wilson's operator product expansion (OPE) [0]. In particular, 
it became possible to quantify the effects of the confining domain on the inclusive 
decay rates. This theory is in the mature stage now (see Refs. |Q, ^ and references 
therein). 

Among the general statements derived for the heavy quark decays, we mention 
here 



• Absence of Aqcj:, /rriQ corrections to all types of fully inclusive decay widths 

• The leading nonperturbative corrections arise in order I/itlq and are given by 
the expectation values /i^, fiQ of the two heavy quark operators, kinetic and chromo- 
magnetic. While the first effect is universal amounting to the correction — /i^/2mQ, 
the Wilson coefficient for the second one depends on the considered process. Both, 
however, are insensitive to the flavor of the spectator(s) ("flavor-independent" cor- 
rections) 0, 1^. 

• The widths are determined by the short-distance running quark masses mQ{fi) 
0]. These are shielded against uncontrollable corrections from the infrared domain 
which would otherwise bring in uncertainty Smq/mQ ~ Aqcd/"^q- 

The effects sensing the spectator flavor per se, emerge at the level l/mg [^, |TU], ^. 



They are conventionally called Weak Annihilation (WA) in mesons. Weak Scattering 
(WS) in baryons and Pauli Interference (PI) in both systems. Their magnitudes are 
given by the expectation values of local four-quark operators.o 

For practical applications we should keep the following in mind (for a recent 
dedicated discussion, see Refs. |]TT], |12|, |[): 

- Good control over the perturbative expansion must be established to address 
power-suppressed effects. 

- The consistent Wilsonian OPE requires introducing the separation of "hard" 
and "soft" scales, with the borderline /i serving as the normalization point in the 
effective theory. 

- One has to allow, in principle, for short- distance (small-coupling regime) effects 
that are not directly expandable in the powers of the strong coupling. 

- Account must be taken of the fact that the OPE power series are only asymp- 
totic |1^, and reconstructing from them the actual Minkowskian observable, generally 



speaking, potentially leaves out the oscillating (sign-alternating) contributions sup- 
pressed, in a certain interval of energies, by only a power of the high momentum 
scale. This is compounded by the fact that in practice one can typically determine 
only the first few terms in the power expansion. 

The last item in the list is behind the phenomenon of violation of local parton- 
hadron duality; in many cases it is among the primary factors potentially limiting 
the accuracy of the theoretical expansion. 

In the actual QCD these technical complications are often interrelated. Therefore, 
it is instructive to investigate the OPE in a simplified setting where these elements 



can be disentangled. As explained in Ref. |TT[, this is achieved in QCD formulated in 



1-1-1 dimensions. Additionally, employing the limit Nc ^ oo one arrives at the exactly 



^The OPE for the inclusive widths, actually, is a priori governed by the energy release rather 
than literally mq [[7|. For simplicity, we do not distinguish between them parametrically unless it 
becomes essential. 

^In the context of the heavy quark expansion, local operators have a more narrow meaning 
denoting the generic operator of the form QOQ, with O being a local operator involving only light 
degrees of freedom. 



solvable 't Hooft model where all the features can be traced explicitly. It is important 
that the 't Hooft model maintains the crucial feature of QCD - quark confinement - 
which is often believed to be tightly related to the violation of local duality. Yet in 
1+1 dimensions confinement appears already in the perturbative expansion. 

The 't Hooft model has often been used as a theoretical laboratory for exploring 
various field-theoretic approaches |]I4|. Most recently the OPE for the inclusive widths 



and the related sum rules in the heavy flavor transitions ||15[ were analytically studied 
in Ref . |]1T[ , where a perfect match between the OPE power expansion and the actual 
asymptotics of the widths was found. The known high-energy asymptotics of the 
spectrum in the model allowed us to determine the violation of local duality in the 
inclusive widths at large ttiq. As expected, it obeyed the general constraints imposed 
by the OPE. Moreover, at least in the framework of this simplified model, the main 
features of duality violation could be inferred from the parton-level analysis itself, 
the working tools of the OPE. The suppression of the duality-violating component 
in Thq was found to be rather strong, with the power of l/rriQ, however, depending 
essentially on the particulars of the considered model and the process. 



Ref. |Tl| focussed on flavor-independent effects. To this end it was assumed that 
the spectator quark q^^ has a flavor different from all quarks in the final state, thus 
ruling out both WA and PI. The OPE analysis of these effects is also straightforward. 
Nevertheless, they may be of independent interest for several reasons. 

First, WA and PI represent power-suppressed and thus purely preasymptotic ef- 
fect. In such a situation one may expect a later onset of duality and more significant 
violations of local duality. Since the above effects are numerically enhanced for actual 
charm and beauty hadrons, studying this question has practical importance. 





Figure 1: a) Quark diagram describing the leading quasifree term in the decay width, 
b) Bare quark diagram for WA. 



Another reason to look more closely at the spectator-dependent effects is related 
to the color-flow considerations usually employed in the context of the large- iVc per- 
spective on QCD, and interpreting the OPE predictions in terms of hadronic states. 
In the case of the quasi-free quark decay width or the WA processes one finds a rather 
straightforward correspondence between the OPE expressions and the hadronic con- 
tributions already in the simplest quark picture where quark allocation over the final 



state hadrons is unambiguous (such a description is expected to hold at N^ — > cxd). 
Let us consider, for example, the free parton decay diagram Fig. la. The ud pair is in 
a colorless state and typically has a large momentum g^ ~ ml flowing through it. It 
is then naturally dual to the contributions from the hadronic resonances in the V-A 
channel (in particular when integrated over g^), much in the same way as in e"'"e~ 
annihilation or hadronic r decays. The c quark together with the spectator antiquark 
produces another string of hadronic excitation. Furthermore, the interaction between 
these two hadronic clusters can naturally be small at large ruf,. WA, Fig. lb, looks 
even simpler in this respect; we will discuss it in detail later on. 



Figure 2: Quark diagram for PI in B meson decays. The weak vertices are broken to 
show the color flow yielding the leading- A''^ contribution. 

The hadronic picture of the processes underlying PI a priori is less obvious. Fig. 2. 
The u quark produced in the decay must be slow to interfere with the valence u. The 
large momentum here flows through the diquark loop (cd) which therefore represents 
the "hard core" of the process. The practical OPE, effectively, prescribes to replace 
the propagation of this diquark by a nearly free di-fermion loop, which amounts to 
evaluating its absorptive part as if the production of the free quarks was considered. 
Basically, no distinction emerges compared to the color-singlet qq' pairs in Figs. 1. 
This may leave one with the feeling of discomfort, for no colored states (in particular, 
with the diquark content) is present in the physical spectrum. In other words, the 
diquark conflguration per se cannot be dual to the mesonic states at any arbitrary 
large momentum transfer. 

Alternatively, one can combine a "hard" quark from the loop in Fig. 2 with the 
slow spectator antiquark to have a color-singlet meson-like conflguration. However, 
such a pair naively is not "hard" : at least in the perturbative partonic picture with 
Psp ~ ""^sp -^ its invariant mass vanishes irrespective of rrif,. While such reasoning is 
clearly of the hand-waving variety, it illustrates nevertheless that interference effects 
are more subtle. 

A more troublesome feature of the interference is also illustrated by the observa- 
tion made in the early 90s by Shifman [|l^. He considered a more general scenario 
with both charged- and neutral-current type interactions, as described by the effective 



weak Lagrangian 
G r 



'^=- r- 



ai(c7^(l-75)6)(d7^(l-75)«) + a2{d-f^{l—f^)h){cY{l-l,)u) 



H.c. 



(1) 

The leading (rather than the power-suppressed spectator-dependent) width was ad- 
dressed. The parton result depends on the color factors ai, a2 in the following way: 
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Q ~ Nc[al + al + ^aifls ) • (2) 



On the other hand, the usual counting rules yield the decay amplitudes into the 
two-meson final states in the form 

M oc \ Nc iai + 17^02 ) for "DgTr"" states 



M oc ^iVc ( aa + -rrCii) for "L)ir" states, (3) 

where, for illustrative purposes, we call gsp the strange quark to simplify distinguishing 
between the two different ways to pair the quarks into mesons. (Since we discuss the 
leading free-parton amplitude, the flavor of the spectator is chosen to be different 
from all other quarks in the process.) Adopting the rules Eq. (0) one gets 

Vq ~ N,l^al + al + —a^a^ (4) 



'c 



more or less independently of the dynamics. While the dependence for the terms 
~ a{ and ~ a^ is reproduced, there is a clear mismatch between Eqs. (H) and (^) in 
the term describing the interference of the two different color amplitudes |T^ . 

There is little doubt that the formal OPE asymptotics must work at arbitrary A^^- 
The arguments above might suggest, however, that the onset of duality is delayed for 
suppressed effects, for example, grow with A^^- 

In reality, we do not think that there is convincing evidence supporting such 
reservations about applying the OPE to flavor-dependent corrections. To provide an 
additional justiflcation, we have explicitly analyzed both PI and WA in the 't Hooft 
model. We have found complete consistency with the OPE, with the onset of duality 
largely independent of the details. As a matter of fact, the parton-deduced OPE 
expression for PI appears to be exact in the chiral limit when all involved quarks 
(but Q) are massless. The resolution of the above paradoxes emerges in a rather 
straightforward manner as well; we will comment on them in subsequent sections. 

We note that we disagree with the claims of the recent paper |]T^ which found 



a mismatch between the actual WA width and the OPE-based prediction, relying 
on numerical computations. We have determined the leading effect analytically and 
showed it to coincide with the OPE result. We comment on the apparent drawbacks 
in the analysis of Ref. Jl^ in Sect. 6. 



The paper is organized as follows. After this introduction, in Sect. 2 we sketch the 
aspects of the 't Hooft model important for addressing weak decays. In Sect. 3. we 
analyze the effects of WA at A^c ~~* oo and analytically compute the large-mQ asymp- 
totics of the corresponding heavy meson weak decay width, with technical details 
given in Appendix 1. Sect. 4 addresses PI; we analytically compute this width up to 
terms like Thq /'ttiq and find full agreement with the expressions obtained in the OPE. 
The effects of local duality violation at large tuq are quantified. The special case - 
with massless final-state quarks - is identified where duality violation is totally absent 
from the spectator-dependent part of the width. In Sect. 5 we present a more detailed 
derivation of the total decay width up to corrections l/wig explicitly accounting for 
nonzero light-quark masses, to demonstrate consistency with the OPE (a detailed de- 
scription of this analysis had been omitted from Ref. W^)■ Sect. 6 comments on the 



analyses which have claimed observing inapplicability of the OPE predictions based 
on numerical computations. Sect. 7 comprises conclusions and overlook and outlines 
our perspective on the problem of OPE and duality violation in the decays of heavy 
flavor hadrons. 

Most technicalities are relegated to Appendices. Appendix 2 collects a number 
of relations useful in constructing analytic I/ttlq expansion in the 't Hooft model 
and summing the exclusive widths. In particular, we give simple expressions for the 
leading terms in the transition amplitudes in Appendix 2.2, perform the differential 
fixed-g^ semileptonic decay width summation up to ^/ttlq corrections in Appendix 2.3, 
prove the OPE prescription for the domain of large g^ and demonstrate the proper 
functional form of the transition amplitudes in Appendix 2.4. The expression for 
the IW functions in terms of the 't Hooft eigenfunctions is quoted in Appendix 3. 
Appendix 4 reports a direct covariant computation of the perturbative radiative cor- 



rections performed while working on paper [11 1; it shows that the result coincides 



with what is obtained by summing exclusive decay channels. 

2 The 't Hooft model and heavy quark decays 

The 't Hooft model, the 1+1 QCD with N^ —>■ oo has been described in many papers 
jTSl , p^ , |2y, |2l[]. The first dedicated studies of heavy quarks in the 't Hooft model date 
back to the early 90s |Q ^. Recent paper |11| specifically addressed heavy quark 
decays and the OPE in this model. Here we only recapitulate some basic features. 
The Lagrangian has the form 

^1+1 = -^ G'^.^G^ + Y^i^S I/) - m,)V^, , iD, = id, + AIT<^ . (5) 

The coupling Qs has dimension of mass. With the above normalization of the gauge 
field, A^ still has dimension of mass as in D = 4. The fermion fields V'(a;), however, 
carry dimension of m}/'^ . 

The OPE analysis is carried out universally for arbitrary number of colors, and so 
far Nc is kept finite. Anticipating the large Nc limit for the final analysis, we define 
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a parameter (3 



^■^^|(iV.-^). (6) 



" c 



that remains finite at A^c — ^ oo. It plays the role of the nonperturbative scale Aqcd- 
Following the actual Standard Model, we choose the weak decay interaction of 
the current-current form. Since in D = 2 the axial current is related to the vector 
one, J^ = e^^J^, we simply consider the VxV interaction: 

G 

jC^eak = — /= {ql^lQ) ii'al^i'b) + H.c. , (7) 

where the dimensionless G is an analogue of the Fermi constant. For semileptonic 
decays ilJa,b are colorless (leptonic) fields. In what follows our main interest lies in 
nonleptonic decays with ipa,b being the quark fields. To make the notations more 
transparent, we adhere to the cases of interest in actual QCD and denote the ip fields 
as u and d quarks, while Q will be a synonym of the b quark, and q called c quark 
(whether we chose m^ ^ Aqcd or consider rrig ^ Aqcd)- The spectator quark q^p 
can be either u or d (for studying WA or PI), or different in flavor from both. 

To address inclusive widths of a heavy flavor hadron Hq one considers the forward 
transition amplitude appearing in the second order in the decay interaction [^: 

Th^ = 2Im |d^X^^(/fQ|zT{£„eak(x)£teak(0)} I^^q) • (8) 

In the limit Nc -^ oo, with Hq being the mesonic (Qqsp) states, factorization of the 
amplitudes holds, which takes the following form for the transition operator: 

|d^x(i7Q|zT{£.eak(x)4eak(0)}|i^Q) = ^ / d^O; T'-(x) n^,(x) , (9) 

where we have introduced the "semileptonic" T^^^ and "hadronic" 11^,^ tensors: 

n^,(x) = {0\tT{d{xh^uix)uiO)-fJiO)}\0) , (10) 

T^'^(x) = {HQ\zT{q{x)rQ{x)Q{0)Yq{0)}\HQ) . (11) 

The Cutkosky rules then yield 

^Hq = ^'j^ fd''xlmT'^''{x)lmU,,{x). (12) 

The factorized representation of the decay width holds only at A"c -^ oo where 
the momenta of the '?/^a'?/'b-pair and {qqsp) become observables separately. In other 
words, in this limit there is a rigid quark allocation over the particular hadronic final 
state and factorization of the corresponding amplitudes, and there is no "cross-talk" 
between them. Yet, Eq. (|T^) represents a certain observable at arbitrary N^ and, as 
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such, enjoys the full rights of being studied regardless of the details of the model. 
In particular, at large energy release it is a short- distance observable and can be 
subjected to an OPE anatomy. In what follows we will discuss this quantity and 
refer to it as the inclusive decay width as motivated by the large- A'^c limit. 

It is worth noting at this point that the qualitative difference between nonleptonic 
and semileptonic inclusive widths disappears for Nc ^ oo. The nonleptonic width is 
given directly in terms of the differential semileptonic distributions ^p^ (though, in 
D = 2 one may have to consider the decays with massive leptons as well). Indeed, 
with rriu = Tn^ as an example, one has (in the momentum representation) 

Ii,M') = -n(g^) (q'S,,-q,qJ) , pig') = --ImUiq') , (13) 

and 

rS, = / dr P(,^) r.(,^) w.th r.,(,^) = ^ ^^ ' <"> 

In D = 2 the correlator of vector currents for massless quarks is known exactly 
and is very simple: 

n(g^) = ^, p{q^) = NJ{q') . (15) 

With nonzero quark masses the spectral density shifts upward, to the mass scale 
~ I3m or m?. A high-energy tail in p also appears ~ Nc{m^ + m'^)/q^. This will be 
quantified in Sect. 3. 

More specific for heavy quark decays is the "semileptonic" part T^^{x), Eq. dTTI) . 
The general color counting rules determine its N^ behavior: 

T^.(x) ~ N'^ . (16) 

Such a leading contribution, however, can arise only with the vacuum as intermediate 
state; all other contributions scale as N^, or even are further suppressed. The vacuum 
intermediate state is possible only when the decay quark q has the same fiavor as ggp- 
This is the effect belonging to WA. Therefore, one has 

(17) 
at g 7^ gsp • (18) 

Since WA is a leading- A'c effect, vacuum factorization saturates T^^ at Nc ^ oo, 
and the effect takes the simplest form. This is the subject of the next section. 

On the other hand, the "usual" non-spectator widths are formally subleading in 
Nc (even though they may yield the dominant contribution to the decay width for a 
particular type of the heavy meson). For such amplitudes the naive factorization does 
not hold, and the explicit expressions take a far less trivial form. In the context of the 
OPE, this emerges as "color-disfavored" structure of the resulting local operators, so 



nl 


- N^, 


pWA 

J- si 


~ Nc 


r„i 


~ N^ , 


Psl 


~ iV° 



that a priori the factorization cannot be appUed to evaluate their expectation values 
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In the limit A'"c ^ oo the spectrum of 1+1 QCD consists of mesonic quark- 
antiquark bound states which are stable under strong interactions. The meson masses 
are given by eigenvalues of the 't Hooft equation 



MlM^) 



m\ — l3^ ml — 13"^ 



X 1 — X 



^„(^) - /3^ /' dy /^ , (19) 



{y-xy 



where mi2 are the bare quark masses of the constituents, and the integral is un- 
derstood in the principal value prescription. The solutions to the equation are the 
light-cone wave functions (p{x), with x G [0, 1] having the meaning of the portion of 
momentum carried by the (first) quark. They are singular at x = and x = 1 where 
their behavior is given by x^° and (1 — x)'^^, respectively, with 70^1 defined by the 
following conditions: 

7r7o _ ml-p^ 7r7i _ ml - (3'^ 



tan7r7o /3^ tan7r7i /3^ 

In full analogy with nonrelativistic quantum mechanics, the eigenfunctions ipn form 
a basis (complete in the physical space): 



dx (friix)ipkix) = Snk , J2 '^riix)ipn{y) = 5{x - y) . (21) 



The weak decay constant of a particular meson is given by 



fn = \— f'dxM^), (22) 

V TT JO 



and the polarization tensor of vector currents (at rriu = rnd) takes the form 

n(g^) = - E -Sw^ ' ^^^') = - E /n ^ {^' - K) ■ (23) 

n H -'"n n 

As mentioned above, at m„ = nid = one has Mq = and /o = JNc/tc, but for all 
excitations fn = 0- 

The transition formfactors between two mesonic states that define the non-annihi- 
lation widths for large N^, are of order N^. Since the weak quark currents Qq are 
formally of order A*"^^, these formfactors are "sub leading" in the same sense as was 
discussed previously and, in general have a more complicated form corresponding to 
the first order correction in the 1/Nc expansion pO| , |2T| . 



Weak Annihilation at Nr ^ oo 



c 



WA in the decays of heavy mesons becomes possible when one of the quarks produced 
in the weak vertex has the same flavor as the spectator antiquark. We assume q = Qsp, 
in our notations. As detailed in the preceding section, in this case there is a single 
contribution to the transition tensor T^^, proportional to Nc and leading to T^q ~ N^- 
This is associated with the vacuum intermediate state, and is given by S 

T^,(x) = ze-^^^^(xo) {Hq\ Qj,q |0)(0| qj^Q \Hq) , (24) 

(with Pq, denoting the momentum of the decaying heavy flavor hadron Hq), so that, 
in the momentum representation, 

ImT^,(g) = ^(27r)^5^(P-g)(i7Q|g7.g|0)(0|g7MQI^Q)- (25) 

This expression is valid in arbitrary dimension for any choice of the weak current - 
in general, one only must replace Q'^^q by the appropriate quark bilinear. Therefore, 

at A^c ^ oo one has 

^X = ^' ^^ ^^q\ ^^/^^ lo)<oi ^^-^ I^q) i"^n^-(^Hj , (26) 

which is illustrated in Fig. 3. It can be traced that the OPE corresponds to the same 
expression if the expectation values of all the higher- dimension four-quark operators 
reduce to their vacuum factorized values (for earlier discussion of WA in a similar 
context, see, Ref. P^)- The latter formally holds, in turn, at N^ — ^ c)o. 




Figure 3: WA correction to the inclusive decay width in the large-A'c limit. Shaded 
loop depicts the exact polarization operator. 

In D = 2 for pseudoscalar Hq one has (0| qjfiQ \Hq) = ifug^fiuPu- For simplicity, 
we will further limit ourselves by the case m„ = rud- Then 

r^^ = ^/^,M^^p(M|J. (27) 



■^We neglect the contribution of another, two-particle state \Hq{P)Hq{P)), also corresponding 
to vacuum factorization, but yielding the u-channcl discontinuity. 
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Strictly speaking, in practical applications of the OPE, U^^, itself is usually likewise 
expanded in l/mg. Also, the deviation of M]j from mg and/or the values of /hq are 
expanded around their asymptotic values at rriQ -^ oo. Therefore, the sensible check 
of duality for practical OPE in WA in the framework of the large- iVc approximation is 
only comparison of the actual behavior of Im n(g^) at large g^ with its OPE expansion 
obtained from the deep Euclidean domain. 

For massless u and d quarks, the exact polarization operator of the vector currents 
is given by Eq. (|1^); the WA width, therefore, vanishes. A non-zero result is obtained 
if one considers a scalar (pseudoscalar) polarization operator, or if m„ or m^^ do not 
vanish. The absorptive part p(g^) is saturated by the comb of narrow resonances 
with heights ~ A''^ and widths ~ l/N^. Therefore, the formal limit N^ ^ oo requires 
an alternative to point-to-point comparison of the actual hadronic probabilities with 
the parton-calculated, or OPE-improved short-distance expansion, even at arbitrary 
large energies. This implies a certain smearing procedure for the actual hadronic 
probabilities. 

Note that, according to Eq. (|26|) the width - however singular it is - always remains 
integrable around the resonances (see also the discussion below, Eqs. (0), (|29| - p3D ). 
By virtue of the dispersion relations the integral of the decay width is expressed via 
the transition amplitude in the complex plane. This amplitude is regular even in the 
formal limit Nc ^ oo when the resonances become infinitely narrow. 

Smearing enters naturally when one considers the 'imaginary' part ^[n(s) — n(s*)] 
at complex s, somewhat away from the physical cut at s > 0. According to a 
dispersion relation it amounts to averaging the physical cross section R{s) with a 
specific weight, 

i [n(.) - n(.')l = i / di j— ^ «(S) . ^ = ^» + '^ ■ P8) 

One can also use different choices of the smearing function having singularities away 
from the physical cut. 

A similar procedure, in principle, is required for the inclusive decays of heavy 
flavors. Strictly speaking, one must introduce the complex variable uj to study the 
analytic properties of the transition amplitude in question |^, |T2|, |1T| : 



A{uj) = j d'^xe -^^^ {Hq\ I T {£.,ak(x)£Lak(0)} \Hq) . (29) 

It can be visualized as the transition amplitude governing the total (weak) cross 
section of the scattering of a fictitious spurion particle S on the heavy quark, 

S{q) + Hq{p) -^ light hadrons , (30) 

or the weak decay width in the process 

Q -^ quarks (leptons) + S . (31) 
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Such processes would appear if the weak decay Lagrangian is modified from, say the 
conventional four-fermion form to the "four-fermion + spurion" interaction, 



For simplicity, it is convenient to assume, as in Eq. (pQ]) that the spurion field does 
not carry spacelike momentum. 

The amplitude A{u}) has the usual analytic properties, and the discontinuity 
across the physical cut at which the point cu = is located, describes the total 
decay width we are interested in. The OPE for the inclusive widths relies on the fact 
that the short- distance expansion of A{uj) runs in l/(u; — Er) and can be applied 
near the physical point u = exactly as in e^e~ annihilation near a positive value 
of s :» ^QCD- i^r denotes energy release.) To the same extent, in principle, a cer- 
tain smearing can be required if the hadronic probabilities still exhibit the resonance 
structure. 

Thus, there is no theoretical peculiarity in the asymptotic applications of the 
OPE for nonleptonic widths. It does not create a conceptual difference to perform 
a short-distance expansion of a single quark Green function (semileptonic widths or 
deep inelastic scattering), the product of two Green functions (e^e~ annihilation) or 
the product of three quark Green functions (the nonleptonic widths). 

Alternatively, smearing in uj can be phrased as smearing over the interval of itlq. 
Indeed, in the heavy quark limit the amplitudes depend on just the combination 
nih — u, 

A{uj,mQ) ~^(0,mQ-t^) (33) 

(there are power corrections to this relation associated with explicit mass effects in 
the initial state). Therefore, in practical terms one can phrase the smearing as an 
averaging over the interval of the heavy quark mass, which may look more transpar- 
ent. 

After this general digression, we now return to specifically WA in two-dimensional 
QCD. It is commonly accepted that, for the two-point current correlators, both at 
A^c — *■ cxD or finite Nc, the properly averaged absorptive hadronic parts asymptotically 
coincide with the leading OPE expression given by the free quark diagram. As was 
mentioned above, for massless quarks this property holds identically for vector and 
axial currents. For the scalar current the asymptotic correspondence in the 't Hooft 
model has been illustrated already in Ref. [|l9l (for a recent discussion and earlier 



references, see Ref. [^). For the WA width, however, we need the m^^j^-suppressed 



effects. The OPE in D = 2 yields at mu = 7nd = m (for arbitrary Nc) 
n(g2) = N, 



1 2m' . ^/l-^ + l■ 



1 ., In 



Q 



q' g^l-^ ^1-^-1 



2TT{0\'muUU+'mddd\0) /m^lng^' 



q' 



2m 1 „ ,„, _ T,,„v -,, o. .^ I Tnn 

~q' 



Piq^) = Nc—, 7=^^ + 27r{0\muuu + mddd\0)6'{q') + O l^ 




12 



where the first term in both equations is just the free quark loop. There is httle reason 
to doubt the OPE for the subleading terms either. Nevertheless, it is instructive to 
give here the direct derivation of the next-to-leading term ~ (m^ + rn?^/q^ in p(g^) 
directly from the 't Hooft equation. 



We follow here the approach of Ref. [|TT| based on sum rules. In the context of the 
Euclidean polarization operator similar considerations ascend to the earliest papers 



on the model, Refs. [T^, ^. To simplify the expressions, we will suppress the explicit 
powers of N^. which enter in a trivial way, and usually will also omit the mass scale 



factor /5, assuming that all energies are measured in units of (3. Then Eqs. ([2^) , (J2 
take the form 



n(g^) = E 



g2 - M2 



P{q' 



E45(^^-M, 



with 



dxv2n(a;) . 



The completeness of eigenstates yields 



E Cn = E / dxdy¥?„(x)^„(y) = 1. 



On the other hand, integrating the 't Hooft equation from to 1 we get 



dx^„(x) = -— 
/o M^ Jo 

Therefore, we get the second sum rule 

n -^0 Jo \ y 

The integral logarithmically diverges at x 
behavior 



dx 



m: 



m' 



X 



m: 



X 



(fn{x) ■ 



(35) 
(36) 

(37) 
(38) 



m:, 



+ 



m: 



I '■Pniy) = I dx I -^ , 

l—yj JO \ X l—xj 

(39) 

and X -^ 1, which corresponds to the 



U{q 



2. 9'r,°° 1 rnl + ml 



In 



q^ 



+ 



(40) 



q^ q* m^ \q / 

given by the free quark diagram, Eq. (|3^) . The divergence of the sum in Eq. (|39|) is 
associated with the high excitations n. Therefore, quantifying the divergence allows 
one to determine the asymptotic behavior of c^. 

To render the sum in Eq. (|39| ) finite we must introduce an ultraviolet regulariza- 
tion. For the logarithmic divergence the exact way is not essential - one is to add a 
hard cutoff factor d^h^—M"^). For analytic computations the Borel-type regularization 
by the factor e~*^"'^ is usually convenient. 

For the regularized sums (we mark them with the superscript A) the completeness 
condition is modified, 

E '^n{.x)ipn{,y) = G{x,y;A), (41) 
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and the Green function G becomes a "finite-width" 5-hke distribution with the width 

1 
A2 



A ~ i- . (42) 



This regularizes the sum in Eq. (^ 

E Mlcl = r\'dq'piq') = f'dxdy i^ + -^] G{x,y;A) 
^.0. h h \ X \-x) 



» 1 I m . m^ 



M2<A2 



= (ml + mfj (inA^ + const) . (43) 

One has, for instance, for the sum over an interval of highly excited states 

f"; qW P{q') = E Mid = {rnl + m^ In ^ + O (^) . (44) 

^^1 A?<M2<Ai ^1 ^^ ^ 

The sum rule (^31) proves that the asymptotics of the smeared lTiiIl{q^) coincides with 
the free quark loop result through terms m? /q'^. It is easy to see that the nontrivial 
corrections in the OPE also emerge only with higher-order terms in 1/q^. Note that 
Eqs. (|3|-H) hold both for light [m <C (3) and heavy (m ^ (3) quarks. However, for 



the asymptotics to start, the condition A ^ rnu^ must be observed. 

Since J^ = e^^J^ and using equation of motion d^J^ = {rriu + majui-y^d, by the 
same token we showed the leading-order duality between the hadronic saturation and 
the partonic expression for the absorptive part of the pseudoscalar current. A direct 
derivation in the same approach is described in Appendix 1. 

As expected, for large M„ one finds the residues c„ ~ rriu^d- Let us note that for 
light quarks c„ are only linear in mg-. since for light quarks (fn{x) ~ a;™"" at x -^ 
(and likewise at x — > 1), the end points of integration in Eq. (^) bring in the l/rriq 
enhancement. 

Combining the sum rule Eq. (|^) with the asymptotics of the 't Hooft eigenvalues 

we obtain 

Again, these asymptotics are valid if "averaged" over an interval of n. 

It must be noted that the explicit constant in Eq. (^) is not important. A more 
detailed derivation of the large-A asymptotics uses the semiclassical expansion of the 
't Hooft wavefunctions. We show in Appendix 1 that the domain of integration where 
X < j2 or y < -^ yields only a finite contribution to the integral in Eq. (|43| ) (and 
likewise in the vicinity of x = 1 or y = 1). At the same time, in the domain x, y ^ ^ 
the approximation G(x, y; A) = 6{x — y) is applicable. 

With the relation for the WA width Eq. (|26D, the comparison of Eq. (|^) and the 
OPE asymptotics Eq. (0) demonstrates that the smeared width in the 't Hooft model 
coincides with the OPE width at least through terms ml^/rriQ. 
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4 Pauli Interference 

In this section we address the effect of interference in the weak decay width of the 
heavy mesons. As explained in the Introduction, it has an independent interest. 
Similar to the partonic free-quark decay width, PI is a 'subleading' l/N^. effect, with 
r^^ ~ Nc^ rather than A*"^. Therefore, the expressions for the amplitudes are not as 
trivial as for WA. Nevertheless, it is not difficult to demonstrate that, again, the 
quark-based OPE predictions coincide with the actual hadronic widths. 

To incorporate PI we must have the flavor of the antiquark produced in the decay 
of virtual W coinciding with the flavor of the spectator; we shall call it u. Moreover, 
the weak decay Lagrangian must contain two different color structures to have PI at 
the same order in Nc as the free partonic width. So, we adopt, for simplicity, 

>Cweak = — /= (ai(c7^6)(rf7^u) + a2{d-f^h){c-f^u)) + H.c. , (46) 

where, again for notational transparency, we identifled Q with h and called q by c. 
In this case the decay width has three terms, 

Thc, = — AT, ( 0^11 + 0^12 + 2aia2ri2) , (47) 

where Fi and r2 are 0(A'°). Clearly, Ti{mb, rric, rriu, rud) = T2{mb, nid, rriu, rric) holds. 

The asymptotics of the non-interference width Fi (F2) for the 't Hooft model was 

calculated in Ref. |TT[] and shown to be given by the OPE one. Now we address the 



analogous question for F12. 

The leading (in mg) contribution to the decay width described by the free parton 
diagram in Fig. la suggests that F12 ~ 1/A"c. For example, for usual V-A interaction 
in D = 4 one would have 

-p parton -*- -pparton -'- -p parton //loN 

J- 12 = T^ J- 1 = — -- i 2 l4ol 

Nc Nc 

(the explicit factor depends on the Lorentz structure of £weak)- Such A'c-subleading 
effects are rather complicated. This suppression, however, is not always present [^ . 
As discussed earlier, invoking the spectator quark through the spectator-dependent 
effects like WA or PI can bring in an A^c-enhancement by effectively eliminating the 
generic l/Nc suppression of the free quark width. As a result, at the price of a power 
suppression in tjiq one can have the A^c-unsuppressed manifestation of the interference 
of the two color amplitudes in £„eak, 

Fi^) ~ O (Ar°) . (49) 

Thus, on the one hand, studying PI allows one to address the interference of the 
color amplitudes in a straightforward way relying on the 1/Nc expansion. On the 
other hand, considering the term ~ 0102 in the decay width in the limit Nc -^ 00 
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automatically singles out the power-suppressed effect of PI. This goes in contrast 
with the usual situation where isolating PI formally requires subtracting the decay 
width of the similar heavy flavor hadron with the spectator(s) having the same mass 
but with the flavor which is sterile in weak interactions. 

The simple quark diagram describing PI is shown in Fig. 2. To leading order it 
generates the operator 



^pi 



—2aia2 



2K 



1- 



m1 + m^ ' 



{h^,l5u) {u-f^-i^h) - 



2mcmd 



7717 



(bu) (ub) + {bi'^^u) (ui'j^b) 



(50) 



with 2mQK having the meaning of the quark spacelike momentum in the final state: 



K 



{mc + mdf 



{nic - nidf 



Ulr 



rra? 



1/2 



It is worth noting that this contribution is not chirally suppressed. Therefore, it is 
meaningful and convenient to consider it in the limit ttIc = rud = 0. 

For B~ mesons having u spectator, the operators in Eq. ( pO|) have the A'^c-favorable 
color structure and, therefore, their expectation values are given by vacuum factor- 
ization: 



2M 



B 



-{B-\{h,^,u){u^,^,b)\B-) 



2M, 



B 



-{B-\{bu){ub) + {bi-i^u) {ui-i^b)\B-) 



^IbMb 
2(mfe + m„) 



(51) 



In particular, at rric = rrid = one gets 

pPi 



r^2 

-2a^a2 — flMB . 



(52) 



We note that P^^ asymptotically approaches a constant when mq —>■ oo. 

It is interesting that there are no l/wig corrections (at small wic,d) to the above 
result. This is a peculiarity of two dimensions where the absorptive part of the 
(di) quark loop in Fig. 2 scales as the momentum to the zeroth power and, thus, 
does not depend on whether one uses pb or Pb as the momentum flowing into it. 
The corrections to the Wilson coefficient as well as other higher-order operators can 
induce only terms suppressed by at least two powers of inverse mass. 

Let us now consider the decays in terms of hadrons. In the absence of WA, the 
leading- A^c final states are pairs of mesons. The partial decay width B -^ D^tt~ takes 
the general form 



C 



kn 



ml\p\ 



al\AkBn\^ + al\AnBk 



2aia2ReAkAlBlBr, 



(53) 
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where A and B schematically denote the "multiperipheral" B ^ k transition ampli- 
tudes and the "pointlike" meson creation amplitudes, respectively: 



Ak ~ {k\J,^B) , 



B„ 



{n\J,\0) 



(54) 



We denote by p the rest-frame momentum of the final state mesons. The PI term is 
then given by the sum 



^pi 



2aia2^Er^A^:^I-^n- 



8m; 



B k,; 



\P\ 



(55) 





a 



Figure 4: Large-A^c decay amplitudes induced by charge- current (a) and neu- 
tral-current (b) terms in the weak decay Lagrangian. 

Both Fi and F2 are saturated by the final states of the type D^ti~ with various 
excitation indices k and n. However, the production mechanism differs: while the 
"charge-current" interaction ~ ai produces vr" by the weak current "pointlike" and 
D^ in a "multiperipheral" way (see Fig. 4a), the situation reverses for the "neutral- 
current" amplitudes proportional to 02, Fig. 4b. These two sources of the final state 
mesons have distinct features for heavy enough Q: the multiperipherally produced 
mesons have the mass squared distributed in the interval from to ~ Aqcd'^q- The 
bulk of the point-like produced mesons have the mass squared a ^ J^gAqcD or m?. 



Ref. |]Tl| demonstrated these OPE-suggested facts explicitly in the 't Hooft model. 

As a result, interference becomes possible only at a small, ~ ^/mq slice of the 
principal decay channels. This qualitatively explains the l/wig power suppression of 
PI which is automatic in the OPE. 

We will now demonstrate the quantitative matching between the OPE-based cal- 
culation and the hadronic saturation of the interference width. To make the proof 
most transparent, we start with the simplest possible case when all final state quarks 



■*for the vector-like current; it would be evenly spread from about m^ to rriQ if the weak vertex 
were scalar. 
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u, d, c are massless. While not affecting tlie OPE analysis, this limit significantly sim- 
plifies the expressions for the individual hadronic amplitudes, as explained in Ref. [ p]T[] . 



In the case at hand, for example, only n = survives for the decay amplitude ~ ai 
(Fig. 4a) and k = for the amplitude ~ 02 (Fig. 4b). The interference then resides in 
the single final state containing the lowest lying massless D^ and vr". Moreover, the 
corresponding transition amplitudes A between two mesons take particularly simple 
form at g^ = in terms of their 't Hooft wavefunctions ||11|| : 

1 /■! 

^M 777l-(^|eM^^i'l^) = -Qz / dx ipk{x)ipB{x) , (56) 

ZlVl£ JO 

where we have recalled that Bq ~ ie^yfQPl^\ and therefore considered only the rele- 
vant light-cone component of the amplitude. Then we have 



PI ^ G'MbN, 



—2aia2 



Air 



dxipsix] 




2 



-2a,a2 ^ flMB (57) 



(we have used the fact that (po{x) = 1, /o = jNc/rr for massless quarks). The minus 
sign emerges since the direction of the vector playing the role of q is opposite for the 
two interfering amplitudes. 

Thus, the OPE asymptotics Eq. ( |5^ ) is exactly reproduced. Apparently, there is 
no violation of local duality at all for PI in the case rriu = nid = rric = 0! This is 
not surprising - in this limit the only threshold in F^^ occurs at zero mass, and the 
OPE series can have the same convergent properties in Minkowskian as in Euclidean 
space. 

With mu,d,c 7^ the interference effects are saturated by several final state pairs of 
mesons, even if the masses are small compared to ttiq. It is still not difficult, though, 
to check that the leading OPE term Eq. (|5^) is reproduced. We keep in mind that 
at nonzero masses the width exhibits the threshold singularities due to the singular 
two-body phase space l/|p | in D = 2. Since it is integrable, the threshold spikes do 
not affect the width smeared over the interval of mass Attiq ~ l/mq. 

The idea of the proof is suggested by the detailed kinematic duality between 
the partonic and hadronic probabilities. The bulk of pointlike-produced mesons 
have masses squared M^ not exceeding (3'^ or m? , while for multiperipherally-created 
mesons k this scale is ~ fimq or rrisprnQ. More precisely |jll|, for the decay rates ~ af 
(Fig. 4a) 

1 777 

T^ E E ^kn oc -^, (58) 

-^ tot k M„>const-mQ '"'Q 

E E r,. « -L. (59) 



tot Mn^mQ Mk>const-mQ v 

Then, calculating the width we can expand around the free quark kinematics M„ 
Mk = 0. In particular, we set 



W\ Mb [ M2 



l + ^T7^ + ... . (60) 



Additionally, we can expand the transition formfactors in amplitudes Ak in g^ 



MM^J -^MO) + §ml^ 



(61) 



and likewise for ^„(M|). In factoring out l/mg in the slope of the amplitude we 
accounted for the fact that it scales as I/tuq in this kinematics. Indeed, the t-channel 
resonances have masses exceeding itiq, and the kinematics (the fractions of the light- 
cone momenta entering computation of the transition amplitudes, see the next section 
and Appendix 2) likewise depend on g^ only as q'^/mQ. 

To obtain F^^ with an accuracy l/rriQ of the free quark width, we actually expand 
the particular two-body decay amplitude only in g^, that is, do not neglect M|- 
dependence for the amplitude ~ ai or M^-dependence for the amplitude which is 
proportional to 02- The expressions for the decay amplitudes at g^ = are very 
simple [plj] : 



M 



kn 




+ a2-MS) 



^« 



kn 



M^l 



(m; 



I -Ml 



{Ml - Mi 



1 rl 

dxipB{x)ipk{x) ■ dyipniy) 

JO 

1 /•! 

^^ dxipB{x)(fnix) ■ dy(fk{y) 
Jo 



(62) 



Then we have 



2\p\M, 



B 



Ml - Ml - Ml . 



Er 



PI, 

kn ' 



G^MbN, 



-2aia 2 ','"" " E/ dxifB{x)fk{x) / dyifB{y)fn{y) / dz(fniz) / dtifkit) 
k,n 47r f-Jo Jo Jo Jo 



-2aia2 



G'MbN, 
A-kMI 



or 



■^pi 




G^ 



dxipsix) 



-2aia2 ^ flMB ■ 



(63) 



(64) 



We extended summation over k and / in Eq. ([63| ) to include all states, since the 
contribution of additional, kinematically forbidden meson pairs is suppressed by high 
powers of l/mq. 

This expression is valid up to the relative ^/rriQ corrections. Indeed, the leftover 
effect of the slope of the transition formfactor FBk is quadratic in mq/rriQ. For 
example, using representation Eq. ( ^8]) we obtain a sum rule which allows one to cast 
it in the form 



E'^Mn-.M 



«(2) 



d_ 



Y.7r^FBk-M^c 



n ^n 



^y'^B{y)vn{y) -Ck 
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= c. ^^Fsu ■ /' dy My) (^ + ^] (65) 

dq^ Jo \ y l-yj 

(and likewise for 5 Ai}^^). The convergence of the integral over y shows that this 
effect is saturated at small n and is of order -^ or ^^, whichever is larger. 

A more accurate consideration reveals that the two amplitudes in Eqs. (|6^ have 
the factors (—1)" and (—1)'^, respectively, and their product, additionally, the factor 
(— 1)^ "Q " '. (The latter is related to the opposite direction of "g*" in the two 
amplitudes and is readily understood since this is a parity-conserving decay Hq — > 
k + n with the meson parities Phq, (~1)" and (—1)^.) Therefore, the sign of F^^ is 
given by the parity of Hq, which is manifest for the OPE result in D = 2, cf. Eq. (|5TD . 

Thus, we see that F^^ agrees with the expression given by the free quark loop just 
to the accuracy suggested by the OPE. 

It is not difficult to estimate the effects of violation of local duality in PI related to 
the thresholds, for small but nonvanishing mg. Since the two-body phase space is sin- 
gular, different ways to gauge its strength will yield different power of its asymptotic 
suppression. Full information is just given by the nature of the threshold singularity, 
the scaling of the corresponding residues and the asymptotic distance between the 
principal thresholds. This would show the contribution to PI of any new decay chan- 
nel, close to the mass where it opens, where the corresponding width is not literally 
given by the OPE. 

It turns out that the magnitude of local duality violation in PI essentially depends 
on the relation between the final state masses. The strongest effect comes from the 
kinematics where one of the mesons belongs to low excitations while another has the 
large mass close to nib. 

The case when ruu = rrid = (but rric > 0) is somewhat special. Here one of the 
interfering decay amplitudes vanishes at the thresholds, and F^^ simply experiences 
a finite jump: 

q2 q9/2 

l^r^ol - — 2ai as ■ const /b mc — ^d^MuQ-Mk) , (66) 

Mk+i - Mfc - ■ 



2M^^ 



Here we used the semiclassical calculations of the transitions to highly excites states 
given in Sect. 3, Eq. ( ^5[ ) and in Ref. [0], Eq. (79). The latter estimate for the Hq -^ 
D^ transition amplitude is valid up to a factor of order one; accepting it at face value 
would yield -\/37r for the constant in Eq. (|66D. 

The referred asymptotics determined the absolute magnitude of the decay ampli- 
tudes relevant for usual decay probabilities, but not their sign which plays a role in 
interference which can be both constructive and destructive. A more careful anal- 
ysis suggests that the relative sign of the two amplitudes alternates for successive 
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thresholds. Therefore, at m„ ~ m^^ <^ (3 and nic <^ mq we have the following ansatz: 
5Vll ~ ^2aia2- const v^/Bm,^ ^^(-l)^'^ (Mh^-tt/^v^) . (67) 

The amplitude of oscillations in PI scales down at least as l/m^q. At large mq 
the threshold widths are much smaller than even the individual principal widths 
saturating V^^ at mc,u,d 7^ 0. 

When niu^d are nonzero, the picture changes essentially in two respects. First, 
neither decay amplitude vanishes at the threshold, since the two-momentum of the 
lighter meson does not vanish: go — ^1 instead of go = 1^1 — ^Hq — M;hr if Q'^ = 
0. Second, the phase space factor l/\p\ becomes now [2Mi{Mhq — Mthi)]'^^"^ vs. 
l/{MHQ—Mth_T:) for Ml = 0. (Ml is the mass of the lighter meson and its momentum 
is called q here. We assume that Mi is much larger than the resonance spacing 
~ 13"^ /mq.) Otherwise, the scaling of the transition amplitudes remains the same. 
Therefore, in this case we have 

(k) 

thr 

M^S ^ m^ + 7r/3v^. 

Strictly speaking, at ruu^ > 13 additional light meson states contribute, and the 
pattern of the threshold spikes in F^^ becomes less even reflecting the superposition 
of a number of similar structures. Additionally, at rriu ~ mc the individual sign- 
alternating behavior becomes more complicated. 

In principle, with all final-state masses not vanishing, there are thresholds corre- 
sponding to decays where both final-state mesons have masses constituting a finite 
fraction oiinq. The threshold amplitudes for such decays, however, are too strongly 
suppressed, since both interfering amplitudes have the chiral and the formfactor sup- 
pression: 

\Mkn\pi oc \ atfc,n~— ^. 



''- 0^ T ^-- --^^ M|- ?H)^ ; , (68) 



The phase space for such decays is also smaller since \p\ ~ ^mq{MHQ — M^^^y/"^. 
Additionally, these thresholds are spaced very closely, at distances scaling as P'^/rriq. 
Therefore, they are subdominant for duality violation. They are related to the sub- 
series of the OPE terms which appear only beyond the tree-level perturbative com- 
putations. 

To conclude this section, let us describe the physical picture which emerges from 
the analysis. In particular, we can see how the interpretation problem mentioned in 
the Introduction is resolved. As expected, the explicit analysis yielded nothing about 
colored diquark correlator, directly. Instead, we observe the duality of differently 

21 



combined quark-antiquark pairs to the hadronic states: one energetic quark (c or d) 
is to be combined with the 'wee' spectator antiquark or slow u produced in the weak 
vertex. It is the pair of quarks picked up from the different final state mesons that 
corresponds to the large invariant mass in the quark diagram. The completeness of 
the hadronic states - or, in other words, the duahty between the parton-level and 
mesonic states - is achieved already for a single fast moving decay quark when it picks 
up a slow spectator. In particular, the 'hardness' of these processes determining the 
applicability of the quasifree approximation, is governed by the energy of the fast 
quark rather than by the invariant mass of the pair. 

There is nothing wrong with considering the colored diquark loop as nearly free. 
Since the overall color is conserved in the perturbative diagrams, in the full graph for 
the meson decay which would include explicitly propagation of the spectator, there is 
always a color mate for any quark in the "partonic" part of the diagram. Moreover, if 
the leading- A'"c contribution is considered, such a color pairing {i.e., which pair must 
be embodied into a meson) is unambiguous. 

Of course, the invariant mass of a single, even fast on-shell quark vanishes. How- 
ever, this does not make the inclusive probability for it to hadronize by picking up 
spectator and forming a meson, a "soft" quantity. For it is not the invariant mass 
but the (rest frame) momentum that determines the hardness. Indeed, the color of 
the initial static heavy quark Q is compensated by the slow spectator. This initial 
distribution of the color field marks the rest frame and makes the hardness param- 
eter for the total probability to look non-invariant if the final state is considered 
perturbatively as a pair of free partons. 

5 Total (spectator-free) width through l/rriQ 

The inclusive decay widths of heavy hadrons in the 't Hooft model in the absence of 
the flavor-dependent spectator effects were considered in detail in Ref. [O. It was 



demonstrated that the analytic summation of the widths for the accessible two-body 
modes reproduces the l/mg expansion of the widths in the OPE, at least through 
the terms high enough in l/mq. In particular, the hadronic width does not have any 
l/rriQ correction which would not be present in the OPE. 

The analysis was performed for arbitrary rric and rngp, but simplified significantly 
when rriu = nid = was set (in the notations of the present paper). Since there is 
httle doubt that the dependence of the hadronic width on m^ and m^ is suppressed by 
at least two powers of l/mq, this simplification cannot affect the conclusion regarding 
possible non-OPE l/mq terms in the width. Nevertheless, we find it instructive to 
describe the direct computation of the terms ~ m\ ^/itlq in the width based on the 
't Hooft eigenstate problem, following the approach of Ref. [rT| and the analysis of 
the previous sections.0 In particular, it illustrates that the case of nonzero masses is 



^The consideration below was elaborated while working on paper pl| , but was not included in 
the final version for the sake of brevity. 
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not any different from 771^ = 171^1 = 0. 

As before, we assume for simplicity tliat rriu = in^, so tliat the uj^d current is 
strictly conserved, and represent the large- A^c nonleptonic decay width as an integral 
of the differential semileptonic width rsi(g^) over g^ weighted with the spectral density 
p(g^), Eq. (|n[). The upper limit of integration q"^^^ comes from vanishing of Tsiiq"^) 
at g2 > (M^ Mo^)2. 



rsi(O) was analytically calculated in Ref. |TT|: 



rsi(o) 



G"^ rriQ — rn^ 



An 



^Q 



mQ 



M. 



dx 



Hq Jo X 



V>H, 



(x) + 0^ 



771'} 



G^ mi) - mi 



An 



n^Q 



,t !M 



(HQlxhmXQlHQ) 



2M, 



Hq 




(69) 



which coincides with its OPE expansion. 

On the other hand, p(g^) does not vanish at g^ > only due to nonzero m^^^. 
Therefore, following the approach of the previous section, we expand p(g^) in m^/g^ 
at large q^ and, simultaneously, rsi(g^) in q^ Itrq at g^ = 0. To this end we write the 
width as 



Hq 



9 

""^'^dg^ p(g^) T^M') = rsi(O) I dg^ p(g^) + /^ dg^ p(g^) (r.^g^) - T^] 

(70) 
Since at g^ ~ ttiq the spectral density p{q^) is explicitly proportional to 7n'^/7nQ, in 
the second integral we can use for rsi(g^) — rsi(O) its leading-order approximation. 
With /dg^ p('?^) = ^c, the first term exactly reproduces the corresponding term in 
the OPE, Eq. (^. 

All transition formfactors are expandable in g^/wig (except, possibly, the point 
g^ —>■ rriq). Therefore, the smeared width rsi(g^) is likewise expandable in g^/m^, 
and so is Thq - at least up to small corrections ~ 1/"^q associated with the domain 
of g^ close to m^. Thus, the second term scales only as rri^^^/TriQ. 

In order to calculate this term, we can use the following facts regarding the 
smeared width: 

• The exact 'semileptonic' width rsi(g^) coincides with the free width r*p''(g^) to 
the leading order in l/mg when mg — -y/g^ ^ p. 

• The smeared rsi(g^) is an analytic function of g^ and is expandable in 1/(??2q — 

• The smeared rsi(g^) does not blow up at g^ — >• tuq. 

We shall comment on them below. Accepting these three facts for now, and neglecting 
m^ compared to tUq we find 
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and 
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Here we have used the sum rules Eqs. (^31) -(^ 
Thus, we get through order l/mg 
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where the kinetic expectation value 
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represents the low-moment um part of the l/niQ expansion of the integral in Eq. 
whereas the term 13"^ /m^ accounts for its "hard" part |ll|]. The OPE result is thus 
reproduced explicitly. 

Neglecting ttIc compared to niq was essential in the above computation. For, at 



rric 7^ the unique linear in mu^ effect appears which was calculated in Ref. [O 
Sect.VI.A : 
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Of course, it has the direct OPE counterpart, see Figs. 5. At rric ~ /? it scales as 

/3^(|mn| + |md|) 





Figure 5: Diagrams responsible for the linear in m.^^^ corrections to ^Hq- 

It is straightforward to obtain the l/mLq term in Thq as well. The least trivial 
corrections not related to ruu , m,d are all incorporated in Eq. (^) . The only remaining 
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part is the l/mq term in the exphcit " ^i '' correction. It is associated with the 

domain of maximal g^, i.e. mq — i/g^ ~ mdS. It has a logarithmic enhancement, 

'^'^ " In ^^ coming from the domain rric <^ rrig — v^ *^ ^^q (and, therefore, 

this log conforms with the one in the free quark phase space). On the other hand, at 
non-equal ttLu and m^ and light c and gsp the smeared width contains the chiral log 
of the form v™M-'|^d) p Jj^ — J — ^ from the domain \ jSirric + ?7^sp) < '^o — V^ < Z^- 
These contributions are given by the corresponding four-quark expectation values 
{Hq\QT ^ccT uQ\Hq) , in agreement with the general proof of Ref. [Q.Q However, 
these operators being A^^-subleading (the expectation values scaling as N^. rather than 
A*"^, unless WA is possible), their values are anyway expressed in a rather ugly way in 
terms of various 't Hooft eigenf unctions both in the s and t channel. Calculating this 
maximal-g^ effect requires straightforward l/rng expansion of the formfactors and 
wavef unctions, which is not too instructive. This is briefly outlined in Appendix 2. 

Now it is time to comment on the assumptions used in deriving the correction 
— "^" '2^'^ . The first one was related to the semileptonic width at nonzero g^. Of 

course, even the stronger statement holds: rsi(g^) coincides with Tl\^'^{q'^) up to the 
terms ~ (3'^ /{uiq — \f(pY. Proving this goes along the same explicit l/rriQ expansion 
of the wavefunctions and the formfactors |Tl|. Namely, to order l/mg one still has the 
transition formfactors determined only by the overlap of the initial ipHg and the final 
ipk wavefunctions. The only difference is that g_ 7^ at g^ 7^ and, as a result, the 
arguments of the wavefunctions change. This purely kinematic modification accounts 
for all changes to order I/tuq. Some technicalities are given in Appendix 2. 

The situation changes at order l/mq. The peculiarity of the point q^ = is that 
the vertices do not renormalize at g^ = 0. More precisely, in the light-cone gauge the 
vertex corrections in Fig. 6a are proportional to g_ and thus vanish, to all orders, 
at g^ = |Tl|. The only surviving contribution is the effect of renormalization of 
the external quark legs. Fig. 6b which is readily summed up to all orders. In the 
light-cone formalism it consists in overall (IR divergent) shift in the reference point 
for the light-cone energy p+, and the dispersion term — /3^/(2p_) formally coinciding 
with the replacement m^ —>■ im? — (3'^ for all quark flavors. 

This perturbative computation has an exact analogy for the actual meson form- 
factors. The absence of the vertex corrections means the absence of the t-channel 
resonance contributions at g_ = 0. This fact is directly seen in the explicit expres- 



^At mu — md these expectation values are multiplied by S^i, — v^v^, with v — Phq/M^q, con- 
taining only the spacelike component, and then the contribution of the lowest pseudoscalar, "pion" 
intermediate state comes proportional to its momentum. From the point of view of hadronic decay 
modes, the pion threshold decay amplitude likewise vanishes at to„ — rud since, for the conserved 
vector current, only the current-meson couplings ~ e^a survive, and the decay amplitudes into the 
mesons with the same parity as Hq vanish at the threshold. These suppressions are eliminated 
when TUu ^ md, and the chiral log appears both in the four-quark expectation value due to the 
pion contribution in the timclikc component, and in the smeared width due to the pion phase space 

~l/(MjjQ-Mthr). 
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Figure 6: Vertex (a) and external quark leg (b) renormalizations of weak decay 
amplitudes, and their hadronic counterparts. 



sions for the formfactors in the kinematics where the fraction of the total light-cone 
momentum a; corresponding to the momentum transfer q, goes to zero. All the strong 
interaction effects occur in the initial and the final states and are described by using 
the exact eigenfunctions instead of the plane waves. 

At g_ 7^ the vertex corrections appear to order g'j, and already the first g'^ 
correction in the coefficient function of the operator QQ becomes different from (1 — 
/S^/mg)"^/^. The vertex corrections are dual to the t-channel resonance contribution 
in the exact amplitude, and they also appear explicitly with the factors g_ and (3'^. 
Additionally, the original tree-level overlap is modified both kinematically and due 
to the additional terms in the expansion of the Q'Juq current in terms of the light- 
cone spinors. All this is directly observed in the explicit expressions for the meson 
transition formfactors. 

Regarding the two other assumptions, the fact that the transition formfactors 
are regular functions of g^ can be seen, of course, from their most general analytic 
expressions. Another way to visualize this is to use the approximate scaling behavior 

{k\Uq)\HQimQ)) ^ -±= {k\J^(q')\HQ(m'Q)) , (76) 
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where q' is adjusted in respect to m'n to have the same rest-frame momentum of the 
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final-state meson k 
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(77) 
This freedom can be used, say, to make momentum transfer q light-like. 
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and then use the exact representation of the amplitudes via the simple wavefunctions 
overlaps. This trick allows one, for example, to see the strong suppression of the 
transition amplitudes to highly excited states M| ~ mg ^ Pttiq at non-zero g^ as 
well, without going into details described in Appendix 2. The relation Eq. (^) is 
violated by the /3^" radiative corrections and by the subleading l/mq terms in the 
expansion of the current J^, which are both power suppressed here. It can be derived 
directly from the explicit solution of the model, see Sects. 2 and 4 of Appendix 2. 



6 Coraments on the literature 



Two recent papers |2^, |T7| claimed to have established inapplicability of the OPE for 
the heavy flavor widths considering the solvable 't Hooft model. Ref. 1^^ addressed 
the conventional spectator-independent decay channels. Numerically evaluating the 
possible two-body decay rates for the values of mg up to mq = 14, it found a small 
systematic excess of the decay width over the free quark diagram which was fitted as 



^Q 
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(79) 



This was regarded as the demonstration of (a priori proclaimed) non-existence of the 
OPE for the nonleptonic widths. Since in the large- A*";, limit the difference between 
the nonleptonic and semileptonic widths disappears, this - if true - would mandate 
the same absence of the OPE for the semileptonic widths as well, an obvious fact 
ignored by the authors. 

In contrast, Ref. [|l^ accomplished the analytic summation of the large-?7iQ width 
in the 't Hooft model, and no deviation from the OPE was found to the high enough 
orders in l/mq (the exact power addressed depended on the values of the final state 
masses). In particular, the absence of the l/mq corrections to the parton result 
for rUu^d ~ /? -C mq was very transparent. Additionally, the correspondence was 
established between the step-by-step quark-gluon based OPE computations and the 
matching contributions in the integrals determining the meson formfactors in terms 
of the 't Hooft wavefunctions. This allows one to compare the OPE computations 
to the hadronic saturation at the intermediate stages rather than only for the final 
result. 
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What could go wrong with the numerical analysis of Ref. [^? We note that the 
apparent "effect" was really small numerically and, in the fiducial range of mq/ (3 ~ 
8 to 10 constituted only about 1%. Moreover, a closer look at the plots of Ref. P? 



shows that the reported discrepancy somewhat increased toward the upper values of 
mq. As a result, a better fit of the numerical points of Ref. P^ would be achieved 



assuming small 0{l/mPq) (not power-suppressed!) corrections to the width, with the 
numerical coefficient ~ (9(10^^). Unfortunately, a priori targeting the l/mq terms, 
the authors did not explore the possibility of alternative interpretations. 

Incidentally, the scale of the claimed discrepancy lies just in the magnitude range 
of the l/mg effects from the OPE which could be a priori expected at mq ~ 10. 
It turns out that the numerical points for T{rnq) in Ref. [^] can be well fitted by 
the leading-mg expression and adjustable (3'^/m'q terms with the coefficient 0{1). 
Adding arbitrary corrections ~ (3^ /m% with the larger coefficient up to 6 -i- 8 allows 
really perfect fits. Since the actual I/ttiq corrections include the expectation values 
of the four-quark operators which are not given by factorization, they are uncertain. 
Their estimates showed that they are typically significantly enhanced compared to 
the naive dimensional estimate ~ j3. 

Thus, the question would remain open before the actual OPE corrections are com- 
puted. This was accomplished in Ref. |]TT|. While the radiative correction enhances 
the width by the factor 1 + j3'^/[2mq) to order l/wig, the kinetic term —^'^/{2rn^) 
tends to suppress it.l3 Therefore, if, as stated in Ref. |^, the plotted values refer to 
the bare mass mq, it seems that the numerical calculation reported by Grinstein and 
Lebed yielded the result exceeding the asymptotic width by an amount ranging from 
a fraction to a per cent. 

The most apparent resolution, in our opinion, is that the analysis of Ref. p^ 
simply does not control the accuracy of the numerical computations of the amplitudes 
at the required level of a few per mill. This is not surprising and roots to the well- 
known problems of numerical computations. The width at large niq is saturated by 
highly excited states whose wavefunctions oscillate fast. This standard problem of 
the numerical computations of the semiclassical transition overlaps is additionally 
plagued by the general complexity of the numerical solutions of the 't Hooft equation 
and the encountered singular integrals. These problems were alluded to in Ref. |B7 . 



The general lesson one draws from this comparison is not new: it is not easy 
to correctly evaluate the actual asymptotic width in too straightforward numerical 
summations of many exclusive widths, were it a computer-simulated approximate 
solution for a theoretical model with - theoretically speaking - potentially unlimited 
accuracy, or the actual experimental data. The approximations made in practical 
applications usually go across preserving the subtle interplay of different effects which 
underlies delicate cancellations shaping the size and scaling of the power suppressed 
nonperturbative effects. 

The recent paper Ref. [|T7| announced even more drastic numerical mismatch of the 



^The value of /ij for quark masses used in Ref. |2^ has been recently evaluated by R. Lebed to 
be about 0.8/3^ (private communication). 
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actual hadronic width vs. the OPE considering WA, even in the leading order at which 
the effect appears. In Sects. 3 and 4, on the contrary, we analytically computed the 
asymptotics of the spectator-dependent widths and found exact correspondence with 



the OPE. As a matter of fact, Ref. ^7j contains self-contradicting assertions: It was 
stated that the (smeared) current-current vacuum correlator (of the type determining, 
say, cr(e"'"e~ -^ hadrons)) is known to obey the OPE. Simultaneously, the difference 
was claimed established for the effect of WA in the context relying on the factorization 
expression Eq. (0). The starting expression - whether correct or not for a particular 
model - thus equates the validity of the OPE for the WA nonleptonic width to its 
applicability for the vacuum current correlator. Fig. 7 shows the ratio of the (smeared) 
values of the absorptive parts of the correlator, 

Imn^^(g2) 



Imn^-(g2) 



according to Ref. [0|. The variance from unity, in general, does seem to be present. 



Imn(s) 




Figure 7: The actual-to-partonic ratio of the absorptive part of the vector current 
correlator at m„ = m^ = 0.56/3, according to Ref. |^, for N^ = 10 (solid line) and 
Nc = 20 (dashed line). Energy scale is in units of (3. Smearing procedure has been 
applied for the hadronic cross section. 

We do not think that this warrants one to become cautious in applying the OPE 
to the polarization operator. There are certain drawbacks in the computations of 

0. 



Ref. 



For unknown reasons the authors presumed that the width Eq 



becomes 



non-integrable for the large- A'"c theory and, therefore, cannot be sensibly smeared.B 
As a result, instead of considering the current correlator in the 't Hooft model itself, a 



®It is obvious that for any regularization consistent with unitarity, the width remains integrable, 
with the integral around the spike independent of the regularization. This is ensured by the disper- 
sion relation; see Sect. 3. 
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rather ad hoc ansatz was adopted which was intended to mimic nonvanishing widths 
of the resonances at finite N^. The simple resonant representation for n(g^), 



C2 



W) = N^i: 2_u2^- ' -I^ n(g^) = nN^ Jl <i ^ {<f - K) (80) 

was replaced by a complicated model where Im 11 (g^) was due to the two-meson states, 
with their production amplitude containing the resonance terms 



q^-Ml + iMnVn 



At large Nc the residues y/N^Cn ~ y/N^ while the resonance decay amplitudes Ainki ~ 
I/a/A^; the widths r„ ~ l/N^. 

As soon as r„ are small enough, such a prescription does not differ from the proper 
spectral density Eq. (RDI) being a concrete functional choice of the (5-distributions, 



r^ 1 



provided each r„ is saturated by the included decay modes: 

Y,\M{n -^kl)\^-^M = 2M„r„. (82) 



kl 



Here ^ki denote the two-body decay phase space. Under these constraints such a 
model - if not true - is at least legitimate. However, two conditions must be observed: 
First, the resonances do not overlap, r„ -C \Mn±i — Mn\. Additionally, the partial 
decay amplitudes and the phase space factors must be practically constant within 
the (total) width of the individual meson. 

Both constraints are satisfied in the formal limit Nc —>■ oo. However, in practice 
the width grows with the mass of the resonance, while the spacing between the 
successive ones decreases. Therefore, for finite Nc adopted in Ref. fl^ (let alone the 
considered case Nc = 1 !) the first condition was not well respected. Regarding 
the second constraint, the problem gets additionally aggravated by the singular two- 
body phase space in D = 2. Therefore, it seems quite probable that the reported 
disagreement is rooted to the inconsistencies of the model for ImH(g^) adopted in 



Ref. fl7 



A detailed look at the plots displayed there incidentally provides support for this 
conjecture. The plots for (the smeared) Iinll{q^) show a rather unphysical shoulder 
at rriQ ^ 10 which coincides with the mass of the resonance having abnormally large 
width (Fig. 9 of Ref. [[l^])- On the other hand, the agreement between the hadronic 
and quark widths is surprisingly good just at lower masses where the resonances are 
more narrow. The point where the curves start to diverge, apparently shifts upward 
with increasing Nc- 
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Finally, a simpler question remains open about the accuracy of numerical compu- 



tations employed in Ref. [T^. Taken at face value, even adopted sampling of mg when 
only 2 to 3 points fall inside a separate sharp resonance peak, seems insufficient to 
evaluate the smeared width reliably. Whether these effects can explain the observed 
discrepancy ~ 20% at mg ~ 10 to 15, remains to be clarifiedB 

It is worth reiterating that the computational difficulties - quite significant in the 
analysis of Ref. fl^ - to a large extent were a hand-made problem. Both c„ and M„ 



are readily computed without cumbersome triple overlaps involving singular integrals, 
and in any case had been determined to construct the authors' model. Computing 
the smeared width directly from Eq. ( pOD would be then quite straightforward. 

7 Conclusions 

We have examined the inclusive decay widths of heavy flavor mesons in the 't Hooft 
model in the context of the heavy quark expansion, paying attention to the spectator- 
dependent effects sensitive to the flavor of the spectator quark. To the order the high- 
energy asymptotics are calculated, there is no deviation from the OPE predictions, 
either for semileptonic, nonleptonic decays or for the e"'"e~ — >■ hadrons-type processes 
- as anticipated. 

We conflrm that there is no difference for the OPE whether a semileptonic or 
nonleptonic width is considered. What matters is only whether the particular ob- 
servable can be represented as the complete discontinuity of the properly constructed 
correlation function over the cut in a suitable "hard" variable. This variable, a;, is 
the same for both semileptonic and nonleptonic decays, Eq. (^) . After that the only 
difference is the number of the quark Green functions to be multiplied and to per- 
form the expansion of the product in the complex plane. The validity of the OPE, of 
course, cannot depend on these technicalities. 

This does not mean, however, that for all types of decays the predictions of 
practical OPE truncated after the flrst few terms, work with equal accuracy at a 



flxed mass mg. On the contrary, considerations of Ref. |^8|] suggest that the effective 



■'hardness" scale can be smaller than literally the energy release. Accordingly, higher 



onset of duality and larger deviations for nonleptonic widths were obtained in Ref. ||T2 
in the instanton-based model. 

Therefore, in our opinion, attempts to check (let alone to disprove) the OPE itself 
in the concrete model are hardly meaningful beyond illustrative purposes. What has 
a potential of providing useful insights, is studying the behavior of the contributions 
violating local duality, so far the least understood theoretically subject. 



^Since the effect in question by itself is l/q^ in the polarization operator, it could be a priori 
conceivable to have subleading 1/g'^ corrections as large as 20%. However, the OPE ensures that 
the corrections are suppressed by at least two powers of q, and thus must be small; they are readily 
calculated. 
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The "practical" OPE yields the width in the power expansion 

l£» =A„ + A + il+... . (83) 

To tuq m^ 

If the series in l/mq were convergent (to the actual ratio), Thq would have been 
an analytic function of tuq above a certain mass mo pointing to the onset of the 
exact local parton-hadron duality. The actual Thq is definitely non-analytic at any 
threshold (whether or not the amplitude vanishes at the threshold). Thus, the 'radius 
of convergence' cannot correspond to the mass smaller than the threshold mass. Since 
in the actual QCD the thresholds exist at arbitrary high energy, the power expansion 
in Eq. (^) can be only asymptotic, with formally zero radius of convergence in l/mq. 
In practice, the true threshold singularities are expected to be strongly suppressed 
at large energies, and the corresponding uncertainties in the OPE series quite small. 
In actual QCD they are expected to be exponentially suppressed eventually, though, 
possibly, starting at larger energies. In the intermediate domain they can decrease 
as a certain power and must oscillate. 



As was illustrated in Ref. |Tl|], the power expansions like Eq. (^) are meaningful 
even beyond the power suppression where the duality-violating oscillations show up. 
In the case of the heavy quark widths where mass mq cannot be varied in experiment, 
the size of the duality- violating component may set the practical bound for calculating 
the widths. Thus, it is important to have an idea about its size. We emphasize 
that one should always include the leading QCD effects to the partonic expressions, 
rather than compare the actual observable with the bare quark result. In the model 



considered in Ref. [|T^], incorporating the power corrections from the practical OPE 
suppressed the apparent deviations by more than an order of magnitude. 

It is worth noting that we identified the case where the exact quark-hadron duality 
in the interference width is saturated on a single final state. It is realized in the chiral 
limit for the final state quarks. It seems to be just an opposite case to the classical 
Small Velocity limit for semileptonic decays noted by Shifman and Voloshin in 1986 
0. While applicable only for large N^ and in D = 2, the present case is peculiar 
in that the duality is not affected by subleading power corrections. Both cases serve 
as a counter-example to the lore often purportingly equating the accuracy of local 
duality to the proliferation of the final state channels in the process in question. 

The performed analysis elucidates how the general duality between the partonic 
and hadronic widths works out its way at large energies. It can be traced that, to the 
leading order in l/mq, the duality is simply the completeness of eigenstates of the 
hadronic Hamiltonian. It is not even required to explicitly solve the 't Hooft equation 
to establish the leading free quark result for the width, but just to know that the 
solutions form a complete basis.tJ The absence of the non-OPE terms already at 
the next to leading order l/mq is a dynamical fact requiring, for example, a proper 



-"^^The Multhopp technique employed in Ref. p^ does not automatically respect orthogonality of 
solutions with truncation and, therefore, in principle can lead to overestimating even the leading-TTig 
coefficient. 
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solution of the bound-state equations. This enters via the average of the moments of 
the invariant hadronic mass in the final state. 

The 't Hooft model gives an explicit example of the duality-violating effects, at 
least of the "minimal", resonance-related nature. As a general feature, we observe 
that 

— the duality between the appropriately averaged (to eliminate the threshold singu- 
larities) hadronic widths and the truncated OPE predictions sets in numerically well 
rather early, after only a couple of the principal thresholds; 

— the local duality works much better in the heavy quark widths, nonleptonic and 
semileptonic than in a{e'^e~ — >■ hadrons). This applies to both the qualitative be- 
havior (generally, a larger density of resonances implying a more narrow minimal 
interval of smearing) and to the strength of the threshold singularities as well as the 
threshold residues. The qualitative discussion of the underlying reason can be found 
inRef. |2§, Sect. 3.5.3. 



Can these lessons be transferred to actual QCD? Unfortunately, there are some 
essential differences between it and the explored 't Hooft model, which must be 
important for the local duality violation. 

The singular l/\p\ two-body phase space in D = 2 strongly enhances the threshold 
singularities, compared to \p\ in D = 4. In the actual QCD enhancement of the non- 
smooth behavior can rather be expected only from single resonances with masses 
close to Mb- While infinitely narrow at Nc ^ oo, they acquire significant width 
for nf/Nc=l which leads to drastic flattening of the resonance- related combs when 
the width becomes comparable to the distance between the successive resonances. 
Additionally, one expects a denser resonance structure, at least asymptotically, in 
the actual QCD in four dimensions (even at Nc — > oo) than the equal in M^ spacing 
in D = 2. All this would lead to suppression of the duality violation. 

Finally, similar types of the condensate corrections lead to a weaker suppression 
in D = 2 due to the smaller dimension of the corresponding operators. This provides 
a larger room for various possible nonperturbative effects. 

These obvious differences would optimistically suggest that the 't Hooft model 
represents, in a sense, an upper bound for violations of local duality in QCD. While 
this is not excluded, such implications must be regarded with caution. 

In fact, the largest duality violating effects in the nonleptonic heavy flavor de- 
cays can be expected from the resonance structures in the combined (qqqq) channel 
in the final state, embedding the quarks belonging to both the "semileptonic" and 
"hadronic" subprocesses which do not decouple completely. Such states are lost in 
the 't Hooft model. Moreover, the limit Nc ^ oo simply erases the difference between 
the nonleptonic and semileptonic decays in this respect. 

Another peculiarity of QCD in two dimensions is absence of the real gluonic de- 
grees of freedom and of the perturbative logarithmic short- distance corrections typical 
for D = 4: (in contrast to the power-like ones in the super-renormalizable D = 2 QCD). 
So far we have no clues if the onset of duality for excitation of the gluonic degrees 
of freedom is the same or noticeably larger than for the processes describing the evo- 
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lution of the "valence" quarks. Correspondingly, while observing a good qualitative 
duality between the bare quark computations and the actual hadronic probabili- 
ties already at energies ~ 1 GeV, one may have to ascend to higher energies in the 
Minkowski domain to reach quantitative agreement at the level of the perturbative 
corrections. Let us note in this respect that the 'real gluons' in jet physics are reliably 
observed in experiment only at rather high energies. Unfortunately, two-dimensional 
theories do not allow one to study this interesting aspect of local gluon-hadron duality 
even in the accessible model settings. 

As emphasized above, we found that the OPE predictions universally hold for the 
inclusive widths. In this respect, the origin of the paradox in mismatch between the 
size of the interference term on the quark and meson languages, Eqs. (0) and @, is 
quite transparent. The mismatch emerges in the A^c-suppressed term. Addressing 
such effects requires a more accurate control of decay amplitudes beyond the leading- 
Ac counting rules. In particular, it is necessary to account for 1/Nc corrections in 
the color-allowed amplitudes rather than use oversimplified prescription like Eq. (j^) 
(a possible resolution was conjectured already in Ref. |[16[] ). 

A closer look reveals that the problem originates just on the quark configurations 
where all four quarks in the intermediate state have the same color. This configuration 
is l/Nc suppressed and can be simply discarded for the leading- A'c probabilities. 
Within this color combination, there is an ambiguity of allocating the quarks over 
the colorless mesons which is absent otherwise. It is not difficult to see that the naive 
computations based on the amplitude prescription Eq. (^ amounts to counting the 
two possibilities independently, which is not justified a priori. Moreover, it is a clear 
double counting at least in the simple quark picture. To phrase it differently, the 
bases of the final states used in the naive color rules like Eq. (|^), are non-orthogonal 
beyond the leading order in Nc- 

With the complicated confinement dynamics, we cannot know beforehand what 
are the actual decay amplitudes to the particular hadronic states in this case - and, 
additionally, the simple two-meson picture of the final state must be, in general, 
extended when going beyond the leading order in Nc. However, it seems obvious that 
the naive prescription have little chance to be true in a complete theory if it violates 
general requirements even in the simplest case of almost free constituents. 

This, incidentally, is an illustration of the fact that the naive factorization of 
decay amplitudes must be violated at some point at the 1/A"c level - either in the 
corrections to the color-allowed amplitudes, or at the leading level for color-suppressed 
amplitudes, or in both. A dedicated study of the above inconsistency was undertaken 
in the framework of the nonrelativistic quark model in Ref. ||3^ . It was found that the 



problem is indeed resolved there in this way, with the concrete dynamical mechanism 
dominating the modification of the 1 /Ac-suppressed amplitudes depending on details 
of the model. 

In Sect. 4, considering the spectator-dependent preasymptotic correction we for- 
mulated the problem of calculating the interference width in the way that allowed us 
to avoid the complications associated with the analysis of the A'c-subleading transi- 
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tion amplitudes, and to study it in the theoretically clean environment. Incorporating 
the spectator quark via PI made interference an effect appearing in the same order 
in Nc as the non-interference widths. Then the hadron-based computation could be 
performed consistently, and the OPE prediction was readily reproduced. 
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Appendices 

Al 't Hooft wavefunctions in the semiclassical regime and 
the UV divergences 

Here we outline a more accurate determination of the scaling with A of the smearing 
parameter A in Eqs. (p])- (|4^) . 

At M„ ^ P, rriu^ the solutions of the 't Hooft equations are nearly free and 
approach the massless solutions |T8|, pl| 



Lpn{x) ~ V2C0S7™X, M^ ^ TT^/J^W . (Al.l) 

The first expression holds outside the end-point domains a; ^ and x ^ 1 bounded 
by the "classical turning points" 

X, 1-x^— ^. (A1.2) 

n 

We again imply the simplified case of equal quark masses m^ = m^ = m. The semi- 
classical wavefunctions ( [Al.lp determine the ultraviolet-singular part of the Green 
functions. 

We shall adopt the exponential regularization of the sums. Then the Green func- 
tion G(x, y; A) in Eq. ( PD has a meaning of the Euclidean Green function for the 
't Hooft light-cone Hamiltonian, with 1/A^ corresponding to the Euclidean (imagi- 
nary) light-cone "time" zx+, and x, y being "coordinates" in the space of the light-cone 
momentum fractions. Formally, the problem is equivalent to usual one- dimensional 
quantum mechanics on the interval [0, 1] with Hamiltonian given by the r.h.s. of 
Eq. (|19]). The first (local) terms play a role of the potential whereas the integral term 
includes the counterpart of the kinetic energy which is now approximately 7r|A;| for 
the "momentum" k much larger than /3. 
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The regularized summation is straightforward and yields literally for G{x, y; A) 



sinhe 



sinhe 



2/Q2 



4 sinh^ I + sin^ f (x - y) 4 sinh^ f + sin^ ^{x + y) 






(A1.3) 



Since the small-n eigenf unctions can be quite different from Eq. ( [Al.l|) , there can be, 
a priori, additional regular (e-independent) terms in Eq. ( |A1.3| ). However, since at 
e ^ the limit of G must yield ^-function, they are absent. 
The first term in Eq. ([A1.3|) clearly yields 6{x — y) aX e —* 0: 



1 



sinhe 



4 sinh^| + sin^f(a;-y) 
with the width A at finite e amounting to 



+ n'^{x — yy 



(A1.4) 



A 



e 

TT 



A2 



(A1.5) 



The second term in Eq. ( A1.3| ) seems to appear when both x and y are small, 
x,y ^eoil — x,l — y ^ e. In this domain, however, wavefunctions <^„(x) with 
n ^ A^ are different (suppressed), and in reality this term must be discarded. In 
practice the domain of such small x must be treated separately, which is illustrated 
below. Therefore, we adopt 



G{x,y;A) 



1 



sinhe 



4sinh^f + sin^f(x-|/) 



1 1 ,, , 

at — <^x,y <^1- — . (A1.6) 



To address the end-point domains x or (1 — x) ~ 1/A^ we need to account for 
the corresponding behavior of (pn{x) in the "classically forbidden" domain, using the 
language of ordinary Quantum Mechanics. In these domains the wavefunctions are 
suppressed: 



V^nix) 



(t) 



71 



const 



l-x 
1—xt 
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X ^ Xt <^l 

1— X<^1— Xt<^l 



Xt 



m\l3^ 



A2 



Now we consider the regularized sum rule Eq. (^31) : 

/^(A^) = r qMq^e-'^"/''"p{q^) = [' dxdy (— 
Jo Jo 

Using G{x,y] A) in the form Eq. ( |A1.6| ) we get 



— +. 
X 1 



m 



X 



(A1.7) 



G{x,y-A). (A1.8) 



/i(A^ 



2-m^ 



, A2 

m —-— + const 

7rp^ 



(A1.9) 
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where the factor 2 reflects the contributions of both x —>■ and x —^ 1. The logarithm 
is saturated at 

(and likewise with x ^>- 1 — x, y ^ 1 —y)- In this domain the expression for G{x, y; A) 
is legitimate. 

Now we show directly that the domain x ^ /?^/A^ yields only a constant, that is, 
a nonsingular contribution at A — ^ oo. First, we note that Eq. ( pH]) ensures an upper 
bound on c„, const /M^. More precisely, 

m^ 1 1 

\Cn\ < const — -r OC - . (Al.ll) 

7 M^ n 

This, of course, follows also from Eqs. (|4^ , ^5]) .til Then we can bound from above the 
small-x contribution as 

/ dxdy 1 (~^[x,y : A) < const — > / dx — [it nx] 

A^-V \ X l-x 'J —:', n Jo x^^ 



'X^ 



^4 oo g-7r2/32n/A2 / j^ \ 7 ^4 



Here we assumed the cutoff in the integral over x at ~ 1/A^, up to an arbitrary 
constant. The dimensionful factor can be either /3'^ or m^, whichever is larger and 
determines the position of the "classical turning" point. The key property is the con- 
vergence of both the integral over small x and the sum over n due to the suppression 
of ipn{x) in the "classically forbidden" domain. 

Thus, the translation rule for the cutoff energy into the regularization Eq. (^21) in 
the integral over x, y is obtained directly from the 't Hooft equation. 

Using the similar technique, it is easy to establish the asymptotic duality directly 
for the pseudoscalar current correlator 

n,(g^) = N^Y: t7^ , ^Imn,(g2) = iV, ^ dlS{M'^-q') , (A1.13) 

where 

<i..= /|(n|*7.«|0)= i|dx(^ + i^)^„W (A1.14) 

1 /"i , im'^ m. 



niu+nid Jo 



r^ I 771 7T7 \ 

/ dx 1 — (fnix) (even n; for odd n). 

JO \ X 1—x 



-■^^While the exact coefRcient in Eqs. (p3|)-(^) is a subtle thing addressed here, the power of n 
itself is simple enough. 
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(The first form is a direct representation, whereas the second expresses the pseu- 



doscalar density as the divergence of the axial current. A certain parity relation [19 



ascending to 't Hooft ensures they both hold.) To find the asymptotics of dn we 
consider the logarithmically divergent sum 



dg 



,2Pp(?^) 



q' 






„ 2{mu+md)M^ Jo \ x 1—x 
Eq. (^8]) allows one to rewrite it as 



^n[x) I dy 
'o 



m' 



m: 



y i-y, 



My)- (A1.15) 



E 



dl 



1 



M2 2(mu+md) 



?^^ 



ma 



X 



771, 



Vn{x) / dy^y) , (A1.16) 

x/ Jo 



which returns us to Eq. (|39D . 

A2 ^/ttlq expansion of the heavy quark weak decay ampli- 
tudes 

As discussed in Sect. 5, the perturbative corrections to the weak decay vertex ap- 
pear to order /?^. In the light-front formalism they vanish at the kinematic point 
g^ = |ll|], since it can be realized as the configuration with g_ = 0, for which 
no physical states in the t channel is possible. The loop corrections proportional 
to the dimensionful coupling jS"^ can be, therefore, inversely proportional to mq or 



E, 



rel 



[mQ 



■m„ 



^f(pY- Since at mq ■ 



■nin 



\fq^ ~ /? the process is 'soft', we do not 



consider here this domain, and do not distinguish between the scales of mq and £^rei- 
It is important that no dynamic gluon degrees of freedom exist in D = 2. Therefore, 
in the 't Hooft model vertex corrections separately are perturbatively infrared finite 
in physical gauges. This was explicitly illustrated in Ref. fill . 



As a result, the actual transition amplitudes {}z\J\iiq) up to terms P"^ /rriq must be 
given by only some overlaps of the initial and final state wavef unctions, at arbitrary 
(though not too close to {mq — mqY) values of q^. We assume that the initial state 
Hq is either a ground state or has a finite (not scaling with mq) excitation number. 
The final state k can be arbitrary. We will show that the explicit expressions for the 
transition amplitudes in the 't Hooft model exhibit this parton-deduced property. 

Let us consider, for example, the representation used in Ref. P^, although we put 
it in a slightly modified form similar to that of Refs. B^, |22|. One introduces the 



kinematic variable uj which depends on q and the final state meson mass: 



00 



1 + 



Ml 



Ml 
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Mi 



Mh 



+ 



Ml 



Ml 



(A2.1) 
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which, in the hght-cone formahsm has a meaning of the fraction of the momentum 
seen in the infinite-momentum frame carried by the particle with mass ^Jq^ in the 
two body decay of B meson, if another particle has mass M^. This fraction has two 
possible values corresponding to the two possible directions of meson k in the rest 
frame. We chose the above branch to have g_ ^ as g^ ^ 0, as in Ref. |T1[]. The 
light-cone fraction uj has a very simple meaning in the rest frame as well: 

{\-^)MnQ = {\Pk\ + Ek)c.m. , ujMhq = {-\q\ + go)c.m. • (A2.2) 

Two cij-dependent Hq D^''^ overlaps are then considered, 

Ck{uj) = / dyipkiy)^HQ{l-il-uj)il-y)), 

u Jo 

Vk{uj) = -uj dy — ^ A2.3 

Jo y 1- {l~uj){l-y) 

which correspond to the effects unrelated to vertex corrections. 

Additionally, in the t channel of the decay process one can have various bound 
states of Be mesons; we reserve index I for them, and their masses will be denoted /x/. 
For each t-channel resonance there is an cu-dependent triple-meson (^^HqBJ^PD^''^") 
overlap 



(A2.4) 
With these notations one has for the transition amplitudes 

Cm^^m (^I qiuQ \Hq{P)) = -y= J2 2 T^ fiFikiuj) - g^Cfe(^) + niQmgVkiuj) , 

(A2.5) 

q, {k\ q^,Q \Hq{P)) = ^ Y: '"^'^^'^l^' m^ico) + q'Ck{u;) + mQm,Vu{u) . 

(A2.6) 
Here, for convenience, we wrote the invariant combinations of the amplitudes instead 
of the two Lorentz components of the current separately. For example, the decay 
amplitude B -^ D^'^^ +7r*^") in the considered case of vector-like interactions takes the 
form 



Y^ "^ 2 -CiFik{uj) + Pnq^Ck{uj) + mqmqVkiuj) 

. I (1 f^i 



Mkn = —?== Cn 

V 2,11 

(A2.7) 
where P„ is the parity (—1)"'+^ of the tt^"^ state. Here and below throughout this 
Appendix we suppressed the factors N^ and put dimensionful /3 = 1. It is easy to 
see that at g^ = (u; ^ 0) only the term ~ Ck survives and this expression indeed 
reduces to Eqs. ([56| , |6^) . 
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We shall show that the first term in this relation associated with the dynamics 
in the t channel, is l/mg-suppressed compared to the terms given by overlaps C and 
V. This implements the approximate "on-shell" condition for the heavy quark Q; in 
this case it can be accomplished employing the l/rriQ expansion of the initial-state 
wavefunction ^Hq- 

A2.1 Nonrelativistic expansion in the 't Hooft equation 

To analyze the heavy quark system, it is advantageous to introduce the nonrelativistic 
variables:c3 



M„ 



mq + er. 



t = (1 — x)mQ and "^nit) 



V^ 



V^n 1 



in terms of which the equation takes the form 



en + 



6^ + 1 



^n(t) 



<p-l 



+ 






21 — ^ 

mq 



^n(t) 



1 



2 Jo (t - s)2 



2mQ J "' ' \ 2t 

The asymptotics of \l'„(t) at 1 <^ t ^ mg is given by 

with Fn = ^jTvm^fn = / dt^!n{t) . 



(A2.8) 



. (A2.9) 



^nW - § 



(A2.10) 



One can extend the bound-state problem Eq. ( [A2.9| ) from [0, rrig] to the whole interval 
[0, oo) if the exact linear potential term | ((1 — 1/''^q)/(1 ~ V'^q)) is replaced by 
its formal expansion in t/rriQ. This is justified literally only up to ^/toq terms {l/m?Q 
in the wavefunction), when the ultraviolet divergences of the static heavy quark 
expansion first show up. Further terms can be obtained using the explicit large-t 
asymptotics in Eq. ( |A2.10D and the end-point behavior of \l'n(t) at t ^ niq. For 
example, the improved nonrelativistic equation takes the form 



e„ + 



e +r 



2ra. 
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^n(t) 



m^„-l , ^ 1 



2t 



t+ 1- 



t' 
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^n(t) 
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ds^„(s) — 



(A2.11) 



Eq. ( |A2.11| ) can be viewed as a usual variational problem for the non-local Hamilto- 
nian defined as 



(*|H|^') 







dt^(t) 



m2 -1 1 / 

-^ + - 1- 

2t 2 I 
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1- 
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t' 



t' 



TTlr 



1717 



2mc 



2mj 



^\t) 



^^This standard for the infinite momentum frame procedure in the context of the 't Hooft model 
was first considered in Rcfs. |E3, 0- 



40 



- i r dt r ds^p^ . (A2.12) 

2 Jo Jo {s -ty ^ ' 

For example, the analogue of the nonrelativistic expansion for the heavy hadron mass 
takes the form 

,. ,,. 1 , S{e)-{ty 4(t3)-3(t)(t^) + (t)3 (t) 
M^,-mQ = (t)-^+ 2^^ + ^^ ^ 

with the average (...) defined in the standard way as the integral over t with the 
weight |\E'(t)p. Here all averages are calculated with the finite-?TiQ ^(i)- To exclude 
the Coulomb interaction term we used certain relations which are derived in the way 
analogous to the virial theorem in Quantum Mechanics. Namely, for any eigenfunc- 
tion \E'„(t) we can consider the average of 7i over the trial function vA^n(At), and 
require a minimum at A = 1 (a similar trick was used in Ref. [^ for the case of 
infinite mq). In this way one obtains 




K - 1)(7) - (^) - f 1 - iV^) - 2 ^ - 3 ^ = ^ f i) (A2.14) 



,2 _ 

t' ^' ' \^ ''^qJ ^ ' ""^Q "^Q V"^^ 

while 

^ e^+l 1 / 2 ^^,V ^^ A , 1 / 1 \ ^^, l{t^) 



2mn 2 



(K-i)4)-(v>)+i(i-^)«+i 




where (V) denotes the expectation value of the integral term in the equation. 

The bound-state-independent term —l/2mQ in Eq. ( |A2.13| ) deserves a special 
note: it represents the short-distance renormalization of the bare mass we used, 
originating from momenta ~ mq. There is no infrared part of the mass: it enters at 
the scale much larger than Aqcd- The effect of smaller momenta is described by the 
nonperturbative 't Hooft wavefunction rather than infrared-divergent perturbative 
diagram. 

Let us note that the operator t/2 is associated with the breaking of scale- invariance 
of the static 't Hooft equation {i.e. \E'(At) is the solution of the equation with the 



t" ji 5 'T-c 



:t— Tit- 



linear term \H/2). In other words, the commutation relation holds 
where Ti^ is the static 't Hooft Hamiltonian, Eq. ( |A2.12| ) with mq = oo. This fact can 
be used to obtain the variation of (t) when including perturbations, in the analogy 
to the case of QCD (Ref. fll. Sect. II): 

5{t) = -l{5Hi) 

for the perturbation dTii which is a homogeneous rank-/ functional of t. This follows 
from the usual operator relations like {iT{\H, A\,B}) = {[A,B]). For instance, to 
order l/mq one has (t) = {t)oo — {t^) /mq, so that in Eq. (|A2.13|) Mhq — rriq ~ 
(^)oo + {t"^ — i"^) / i'^fng) , as it should be. 
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A2.2 Nonrelativistic expansion of the decay amplitude 

In the l/mQ expansion of the Hq -^ D^''^ amphtude we pass from (fHgix) to the 
nonrelativistic wavefunctions "^Hqit) and, in the case of the triple vertices Fik, rewrite 
if I "(x) via the corresponding \E'; "(u). The nonrelativistic approximation for c quark 
is not employed here, however. Then 

Cfc(u;) = -^ / dt ^nS) V'M 1 - J. ^ , (A2.16) 



iTr I -i-\ ^ 



P„M - -^-^r^ f ' " ^4^w«„m "\^ r^- . (A2,17) 



{l-uj)mQ 



{l-uj),/m^ Jo l-t/rriQ 

and, for the triple vertex, 

Fik{ijj) = uj{l-uj) X 

Hdu /'dy-^^£M^^{vI/^,(^«+(l-^)mQ)-vI/^,((l-^)(l-y)mQ)}. 

(A2.18) 
At arbitrary g^ we have a; ~ 1; nevertheless, at large enough energy release (1 — 
uj)mQ ^ 1 still holds. This parameter defines the 'hardness' and is used in the l/mg 
expansion. 

The nonrelativistic expansion amounts to assuming that the support of the ^ 
functions is limited to a finite interval of the argument of order 1, and extending the 
integration over t and u to infinity. With the fiducial domain of integration t, m ~ 1, it 

1 /o 

is readily seen that Ck and V^ scale like mg and lead to the properly mg-behaved 
transition amplitudes. On the other hand, Fi^ ~ ^/fnq and are accompanied by 

— 1/2 

the factors ci ~ tuq . The sum over / is effectively cut off above / ~ mg where /i^ 
exceeds niQ. Altogether, the terms with Fik yield corrections to the decay amplitudes 
suppressed by at least l/mg. 

The leading heavy-quark transition amplitudes governed by the overlap factors C 
and V, Eqs. ( |A2.16| )- (|A2.17|) , exhibit explicitly the proper functional dependence on 
the combination of tuq and q^: the inner product of wavefunctions depends only on 
{l—uj)mQ whose value just fixes the energy (or momentum) of the final state hadron 
in the rest frame, Eq. ( [A2.2| ). In reality, this property of the leading-mg transition 



amplitudes is more general and is not related to the smallness of the perturbative 
corrections. It holds at small energy release as well, where the vertex corrections are 
not power suppressed; this is addressed below in the end of Appendix 2. 

A2.3 "Semileptonic" width Tsi{q^) at arbitrary g^ with l/mg accuracy 

Here we illustrate how the parton expression for rsi(g^) is reproduced at nonzero q^. 
We assume that the width averaged over an interval of rriQ or g^ is considered, so 
that the threshold factors become nonsingular. The analysis differs from the case 
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of g^ = anatomised in Ref. [|r^ only in technical details. In particular, the width 
is still saturated by the states with M| ^ (1 — g^/mQ)mQ, and the summation over 
the final states can be extended to infinity. The decay amplitudes are given by the 
wavefunction overlaps to this accuracy. Let us assume for simplicity that rric <^ rriQ. 
Then only the term ~ Ck survives: 



nq' 



E 



q- 



t 



(l-(^fc)mQ /•(l-(^fe)mQ / 

dt / ds^^(t)^Hc3(s)¥^Jl-- 

Jo ^ V y {l-uJk)mQ^ 



V^fc 1- 



{l-uJk)mQj 
(A2.19) 
The leading term in the expansion of the width emerges if we neglect Mk compared 
to (1 — ct;)mQ. Then Uk = q^/Mfj and 2\f)k\ = {l—u!k)MHQ are /c-independent, and 
we get up to the power corrections 



nq') = ^ 



G^ 1 



An 1—uj Jo 
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dtdsvl>^^(t)vl/^^(s)(5 -— 



Xl-uj)mQ^ 



G^niQ 
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1 + o{ 

(A2.20) 
Account for the l/mg effects requires expanding |pfc| and Uk in Ml/irriQ — q"^), and 
leads to the 'second' sum rule for the average of M| |]ll| . Here is how it works. 
The partial width Tk{q'^) to l/mq accuracy takes the form 



Tkiq') 
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where ujq = q'^/MJj . 
Using the relation 

E Mk fk{x)fk{y) 
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(A2.21) 



(A2.22) 



{x - yf 



and equation ( A2.11 ) to the leading order in l/mg, we obtain 

t 



Y^Ml J ds^H,{s)^k(^l--^ 
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(A2.23) 
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(e = Mhq — rriQ). Hence, we arrive at 
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r(g^) = -MH, 1 + — 
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jdt{t-2e)^l^{t) + 
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m{2,-t)^'H^{t)^H^{t) 




To evaluate the last integral we note that for any f{t) 
Therefore, Eq. ([A2.24| ) takes the form 



(A2.24) 
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where we recalled that e — {t) = O (^) • 

The next, I/ttiq corrections, on the other hand, are not only limited to the term 

2 

~2r^' ^^^ include also those coming from the perturbative vertex renormalization 

and from the sensitivity of the decay width Q —>■ q + (p with m'^ 7^ to the short- 
distance quark mass renormalization m^ -^ rn? — /3^. 



A2.4 Decays at maximal q^ 

At q^ close to the energy release, (mg — niq)'^ — q^ ^ P^tT'Q (or at (mg — nig) ^ /3) 
the decay processes are not hard but proceed over the time intervals > 1//5. As such, 
they in general are sensitive to the t-channel evolution as well. OPE does not allow 
to compute these widths in the short-distance expansion directly even to the leading 
order in l/mq. However, it relates the overall width associated with this domain of q^ 
(i.e., integrated over q^) to the expectation value of the local four-fermion operator: 



pendpoint ^ ^^(^q) " (^M^. 



^,^Ml^ 



2M. 



Hq 



-{HQ{P)\Q^,qq^,Q\HQ{P)) 



(A2.27) 

(this fact was proven in [^ although was used to calculate certain preasymptotic 
corrections to the heavy quark widths since the original papers 0). Complicated 
strong interaction dynamics shows up here as a nontrivial expectation value of this 
operator (for light g; it is perturbatively calculable when q is heavy). In the 't Hooft 
model it essentially depends on the details of the lowest t-channel states. As has 
been already discussed a few times, a smeared decay width is assumed here, which is 
reflected in Eq. ( [A2.27| ): it incorporates the spectral density pirriQ) averaged over an 
interval. With a continuous p{q^) there would be no resonance structure in T{mQ). 
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Referring to paper |^ for the formal derivation of Eq. ( |A2.27| ), here we illustrate 
it schematically in a transparent way. Let us consider a particular meson k in the 
final state. The corresponding partial decay probability (in the rest frame) is given 
by 

{HQ{P)\Qj,q{0) \k{pk)){k{pk)\ qiMO) \Hq{P)) . (A2.28) 

Integrating over some interval of q^ we can pass to the variable |pfc| according to 



Ml, + q^ - Ml 
2M. 



^0 = WT7 ' ^^0 = -d^fe Ek dEk = \pk\ d\pk 



and have 
Jdq'Mq') = G' j^^^p{q\\pu\))({P-Pk),{P-Pk)u-^,AP-Pkf) X 

-1-{Hq{P)\ Q^.qiO) \kip,)){kip,)\ g7.Q(0) \HQiP)) . (A2.29) 

2Mhq 

The factor of 2 above corresponds to two possible directions of pk and is the unit 
'sphere' surface area in one space dimension. If the smeared width is considered, the 
corresponding smeared spectral density p(g^) replaces p(g^) in the above equation, 
since both the decay amplitudes and phase space depend on the combination mg — y^ 
(for simplicity we imply that g^ ^ ^)- 

With the smooth p{q'^), we can neglect pk compared to Mhq in Eq. (|A2.29|) and 
have, up to l/mg corrections. 



r G"^ 

J dq^ Tk{q^) = ——mlp{nil) X 



{v,v, - V) / 7^^{Hq{P)\ Q^,q \k^){k,^\ q^^Q \Hq{P)) . (A2.30) 

The corrections to this expression appear due to the dependence of p(g^) on g^ near 
g^ ~ rriQ and are power-suppressed. 

If we formally summed over the final states k, the r.h.s. of the above equation 
would become 

-^m|p(m|) iv,v,-5,,) ^^ J ^^{HQiP)\Q^,q{0)\k^){k,^\q^MO)\HQ{P)) 

Hq i^ k 

= ^mlpiml) (v,v,-5,,)^^{HQiP)\{Q^,q)iq^,Qm\HQ{P)), (A2.31) 

which shows the representation of the matrix element as a sum over intermediate 
states. The integration over g^ replaces dq"^ / {27iMHQ\pk\) in the partial decay width 
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by d'Pk/{2TT2Ek) which is just the quantum mechanical summation over the interme- 
diate states. 

Of course, this correspondence is not accidental for D = 2. Indeed, the general 
two-body phase space is given by 

j;A2.32) 
To obtain the decay width we multiply it by the transition amplitudes (HqlQ'j^qlk), 
{k\q'y^Q\HQ) , by :ilmn^^(g) and by {27i)^ /{2Mhq). Integrating over q'^, we get 



Jdq'T.iq') = ^-^ j dq'^2{Ml^,q\Ml) {HQ\Q^,q\k^){k^\q^MHQ) 



= ^' / (o^d'Toj, -l^^AP-Pk),^{HQ\QlAkp)(kp\qiMHQ) ■ (A2.33) 

With Im Yi^u{P — Pk) approximated by its (smeared) value at g = rnqv, this is 
the contribution of the given state k with mass Mk (whether meson, baryon, or 
a multi-particle state) to the expectation value of the local four-quark operator 
{HQ{P)\{Q^,q){q^,Q)\HQ{P)) 0. 

Of course, in reality we do not want to extend the summation to all states with 
arbitrary large energies, which would correspond to integrating over g^ down to q^ = 
— oo. The large-g^ domain is cut at some q^ ^ (ttiq — A)^ ~ rriQ — 2mQA, with 



A > Aqcd- This means that the heavy quark expectation value in Eq. ( |A2.27| ) 
corresponds to the normalization point /i = A. We note that the normalization point 
enters as the cutoff over the energy of the intermediate states, the common case for 
the heavy flavor systems (see reviews 0, |2^ and references therein). 



It is worth reiterating an important point regarding the end-point domain: the 
matrix elements in Eqs. ( [A2.27| ), ( [A2.31| ), (|A2.33| ) are those in the effective theory 



rather than in full QCD - even if there were no high-momentum gluons in the latter. 
The representation of the four-fermion expectation values 

{HQiP)\iQ^,q) iq^,Q)\HQiP)) = E / ^J^D-i2E, ^HQ\Q^,q\k^){k,^\q^,Q\HQ) 

(A2.34) 
in the full theory includes not only the intermediate states without heavy quark Q, 
Fig. 8a, but also those containing the pair of Q and Q, Fig. 8b. The latter are not 
related to the (perturbative) heavy quark loops and appear even in the free theory 
- for example, the state Hq{P) + {Q + q) with Q and q having small spacelike 
momenta ^ A totalling to pk- It can be simply Hq{P) +Hq{p). These states appear 
since the Lorentz-covariant currents QTq, together with operators creating Q contain 
also operators annihilating antiquark Q, which is not the case in the nonrelativistic 
field theory. Such amplitudes are shaped by the time scales ~ 1/A and are as 
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nonperturbative as transitions into the usual light hadrons. Imposing the cut on the 
virtuality itself does not help to eliminate such spurious states with "valence" Q. 
In particular, integrating out the high-momentum modes of quark and gluon fields, 
or imposing a cut on the gluon momentum (virtuality) in the full theory - while 
generating an effective low-energy theory - would not yield relevant operators, at 
least in this context. 



^ 



-#- 



X 



Hq Hq Hq 




Figure 8: "s-channel" with nq = (a) and "u-channel" with uq = 2 (b) intermediate 
states saturating the four-quark expectation values in full QCD. Both types are de- 
scribed by 'soft' transition amplitudes and contribute simultaneously to the forward 
0{G'jp) transition amplitude near maximal g^. Only the targeted "s-channel" states 
appear in the nonrelativistic effective theory. The two points with weak currents are 
shown separated in space (x 7^ 0) for clarity. 

The extra contributions of Fig. 8b are counterpart of the intermediate states in 
the transition operator which are nearly on-shell for the reverse ordered product of 
currents {q'~ivQ){Q'~i^q) in their time-ordered product, 

%M) = J d'^xe-"^^ zT {q^MO)QlM^)} (A2.35) 

at Xq < 0. They are responsible for the w-channel cut of the Lorentz-covariant 
forward transition amplitude T^;^ in Eq. (pH]). In the kinematics g^ -^ 771%, qq -^ niq 
the transitions to both s channel states with ng = rig = and u channel states with 
nq = Uq = 1 are equally 'soft' and nonperturbative, and cannot be disentangled in 
the single amplitude T^;^(g). 

On the contrary, only the necessary light intermediate states are present when 
the proper nonrelativistic effective theory of heavy quarks is considered, where Q{x) 
do not include antiquark operators, either creation or annihilation. The expectation 
values determining the corrections to the widths must be understood only in this 
sense. 

We conclude this Appendix by noting that the leading-rrig transition amplitudes 
Hq — > k obey the stated scaling in respect to niq and g^, i.e., for the particular final 
state k depend on one combination {m^ + M| — q^)/mQ having the meaning of the 
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rest-frame energy Ek. 

-^{k\qrQ\HQ)^2^^^ ^ -^(fc|grQ'|/7Q,),,2,^^, (A2.36) 

if 

ml + Ml -q^ _ mQi" + M| - g/^ 

vdQ ruQ, 

This functional dependence ensures, for example, that the four-quark expectation 
values discussed above have a finite mg-independent limit as mq -^ oo. 

In terms of the light-cone parameters Eq. ( A2.36 ) says that the formfactors, up to 



certain powers of Jtoq must be functions of (1— u;)mQ, see relation ( [A2.2| ). The ex- 
plicit expressions (|A2.3|) - (|A2.7|) do exhibit this property, although in slightly different 



ways in different kinematic domains. To locate it, we once again can be guided by the 
duality between the various pieces in the hadronic expressions based on the 't Hooft 
equation, and the contributions coming from the corresponding Feynman diagrams. 
This duality exists in the physical gauge and is transparent when the light-cone for- 
malism is used. Additionally, the scaling behavior of various perturbative effects can 
be anticipated beforehand if general OPE facts are considered. For example, the high 
infrared stability of the inclusive widths, together with the fact that no gluon degrees 
of freedom (and, therefore, no bremsstrahlung) exist in D = 2, leads to the finiteness 
of the vertex corrections and allows to estimate their magnitude. 

It has been mentioned already that Eq. ( |A2.36| ) holds at large energy E}^ ^ /5: 
then the amplitudes are given only by the simple overlaps C^ and "Pfc. The integrals in 
Eqs. (|A2.16|) and (|A2.17|) manifestly depend only on (\—iS)mQ if "^Uq^^ are nonzero 
only for t ;^ /3. 

A different situation takes place at maximal g^, for example, for light final state 
quarks. The triple-overlap terms proportional to Fw^ are then of the same order. 
However, in this kinematics Fi^ show the same dependence on {\—LS)mQ. Indeed, in 
this case tu — *> 1, and they take the form 



(A2.37) 
The position of the excited states in the t channel is also driven by the nonrelativistic 
expression [ii ~ mg -|- q with mQ-independent e^. Therefore, the pole factors in front 
of Fik are simply proportional to I/uiq: 

4^ ^ ,,T1-, for q' = (mg - A)^ with A « mg . (A2.38) 

On the other hand, A differs from (l—uj)mQ by just a (Mfc-dependent) constant. 



Efc ~ M^^ - mg + A = eo + A , Ek + ^ El - Ml = (l-cu)mg 
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— 1/2 

Thus, together with the factors q oc mq this term possesses the required depen- 
dence on {l — Lj)mQ and overall power scaling in mq. 

This reasoning may seem inapplicable when both Q and q are heavy quarks, in 
particular, at small uiq — irig. The t-channel B^P wavefunctions are then nonrela- 
tivistic in respect to both quarks, and must be treated accordingly. This, however, is 
a transition between two heavy quarks, and even the classic case of the small velocity 
1^. As mentioned above, in D = 2 the amplitudes are finite even at w 7^ and 
the perturbative corrections are inversely proportional to the heavy quark mass (but 
not to powers of 1/g^ which blows up at the safe zero recoil point). And indeed, it 
is easy to find that Fik in Eq. ( [A2.18| ) are universally suppressed by the mass scale 
Ek ~ {l — uj)mQ whether mg is small or large. 

Therefore, the functional relation Eq. ( |A2.36| ) holds in the 't Hooft model for any 
decay kinematics. 

A3 The IW function and the heavy quark distribution func- 
tions 

The IW function ^ is most simply defined as the flavor-diagonal vector formfactor 

{Hq{p')\Qi,Q\Hq{p))=^(^] {p + p% 

in the large- ?tiq limit. The expression for it quoted below can be obtained using 
the l/mg expansion of the amplitudes in the previous section, if we employ the 
nonrelativistic expansion of the final state as well. Alternatively, we can use directly 
the simple universal expressions Eq. (|56|) at g^ = relying on the scaling relations 
([f6|)-([78|) by adjusting the final-state quark mass: 

~2l 



Mhq, = zMhq with ti; = (ot') = — - — , z = w± Vw^ - 1 . (A3.1) 



The perturbative corrections are finite and vanish at rriQ ^ 00 in two dimensions. In 
this way we obtain 



1+W + \/w^ — 1 JO 



The last expression was used to make representation more symmetric. Here '^ Hqit) 
is the solution of the static 't Hooft equation 
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Expression ( |A3.2| ) is easily generalized for the case of inelastic transitions: 

«.-(-) = ^ r d' < (75) < (V5<) . (A3.4) 



This agrees with the expression obtained in Ref. |2^. We found a direct proof (it will 
be presented elsewhere) that the right-hand side of Eq. ( [A3.4|) does not change (up 



to a sign) if 2; ^ 1/z. This ensures the same result for both solutions for z in terms 
of w, the last of Eqs. (|A3.1|) . 



It is interesting to note that the absolute maximum for the (elastic) IW function 
which would correspond to the structureless point-like heavy flavor hadron, is given 
by just the factor EJ 

U..H = \^^ = ^ . (A3.5) 

V 1 + w I + z 

It would be saturated by the light cone wavefunction \I'(t) ~ cj \ft. Indeed, one has 

i{z) = \Y^] for ^(t) =t"e-'^*. (A3.6) 

Of course, such a situation is not realized for any actual hadronic state in the 't Hooft 
model. Constructing such a hard-core wavefunction requires a coherent superposition 
of all the excited eigenstates. 

Studying the decay distribution for the weak transitions into light quarks, m^ ^ 
ttiqAqcd allows one to determine the light-cone distribution function. In particular, 
it emerges directly in the decays at g^ = as the differential decay probability vs. go 
or M,\,, g: 

-2 

hadr 



1 drsl(O) mp^oo ,^^ , ^/ Ml 



mQ{Mu^ - mg) F 1 —^ , (A3.7) 



r,i(0)dM2,dr ^^ " ^' V mQ{Mu^-mQ) 



/Vf2 _ yr2 

go = ^-^ '-'''' at g^ 



where rsi(O) as a function of M^^^jj. is 



E rir(o) 



K<K.i. 



Using the leading term in the l/mq expansion of the decay amplitudes, one can 
check that F{t) coincides with \E'|^ (t). To obtain the distribution over M|, it is 



^•^This function is not analytic and rather has a branch point near w = —1 (or q^ = 4toq). This 
is just a reflection of the Fermi statistics of heavy quarks confined in the bosonic state. This means 
that it is impossible to construct a heavy pointlike meson from fermionic constituents, with the 
radius much smaller than the mass. For light constituents this is possible as exemplified by the 
chiral pion. 
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again convenient to use the explicit cutoff-dependence of tlie ultraviolet-regularized 
Green function G{x, y; A), as it was done in Sect. 3. Earlier discussion for the 't Hooft 
model can be found in Refs. [^, P5| . 

The small-x behavior of the heavy-quark distribution function for decays into light 
quarks has, in general, a non-integer power depending on the spectator mass. 



F{x) oc X 



27sp 



tan7r7sp 



mtp - P' 



P' 



(A3.8) 



Just such a behavior was conjectured for the QCD light-cone distribution function in 
the model suggested in [Q. In D = 4, however, the distribution function itself and 
the exponent 27, in particular depend on the renormalization point. 



A4 Perturbative corrections and IR regularization 

The perturbative corrections to the weak decay vertex explicitly depend not only 
on the gauge but also on the employed infrared regularization, even at g^ = (see. 



e.g. Ref. [0, Sect. Ill B). This subtlety is aggravated by severe infrared divergences 
in the individual diagrams in D = 2, which appear when unphysical degrees of 
freedom are introduced (say, in the covariant gauges). While OPE - applicable 
for inclusive widths - ensures that the widths are rather insensitive to the infrared 
contributions, in practical computations it requires carrying out one and the same 
regularization procedure consistently through all stages of computations. A particular 
infrared regularization, on the one hand, is built into the 't Hooft equation in the 
form of the principal value prescription of the Coulomb exchange integral and the 
self-energy terms. On the other hand, the same infrared regularization was employed 
in establishing the short- distance non-renormalization theorem for the weak vertex. 
Therefore, it is instructive to obtain the same perturbative correction /3^/(2mQ) 
to the decay width starting with the usual Feynman graphs in the covariant gauge 
routinely used in four dimensions. To this end we independently calculated the one- 
loop perturbative corrections to the quark decay width Q —^ q + (p at m^ = 0, with 

the decay interaction 

G 
>Cweak = -7^qQ(t) + H.c. . (A4.1) 

ZTT 

\n. D = 2 this is a full analogue of the four-fermion (semileptonic) decay with massless 
leptons for vector-like weak currents. 

Since all the individual diagrams are too infrared divergent, the computations are 
simple only in dimensional regularization. All corrections are nontrivial, including 
virtual vertex corrections, mass and wavefunction renormalization and the "real" 
gluon emission width. The sum of all contributions takes the following form (at 
rUq = 0) in arbitrary dimension: 



^pl-loop 
ptree 



Cf 



9l 



(47r: 



m, 



,D-4 

Q 



r 2 



D 

~2 



D-3 (D-l)(D-4) 



D-A 



D-3 
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V^{D/2) 3D^ - 28L>3 ^ 109D2 _ 208D + 160 1 
^ r(|D-2) {D-l){D-^){D-Af /• ^^^-'^ 

It is easy to see that the poles at D = 2 and D = 3 all cancel out (at D = 3 heavy 
quark masses acquire infrared logarithmic divergence as well). The remaining pole 
at D = 4 is given by the ultraviolet renormalization of the Yukawa coupling. 
At D = 2 the above expression reduces to 

^pi-ioop _ Cf^^ r*^*^" , (A4.3) 

which coincides exactly with the result obtained in the light-front formalism and by 
the summation of the exclusive hadronic widths in the 't Hooft model. 

We note that here one has an instructive example: the perturbative calculation of 
the decay width is IR safe by itself, and the first-order correction is meant to describe 
the actual decay width with the l/wig accuracy. Higher orders of the perturbative 
expansion yield l/mg and smaller terms, and the uncertainties in summing the per- 
turbative series are exponentially small in mq. However, this is not the complete 
answer already to the order 1/m^. This illustrates incompleteness of the purely per- 
turbative expansion, even in that the actual nonperturbative effects are not signalled 
here by generic intrinsic divergences of the perturbative series. 
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